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1.  Introduction 


1.1  Background 

In  the  area  of  seismic  monitoring  of  underground  nuclear  explosion 
testing,  the  focus  of  research  has  shifted  over  the  years  depending  to  some 
degree  on  the  terms  of  existing  treaties  and  changing  political  interests. 
Throughout  the  era  of  the  Threshold  Test  Ban  Treaty  (TTBT),  the  greatest 
attention  has  been  given  to  determining  the  size  and  character  of  rather  large 
nuclear  explosion  tests  for  specific  source  areas.  In  recent  years  increased 
concern  over  reaching  agreement  on  a  Comprehensive  Test  Ban  Treaty  (CTBT) 
and  in  curbing  proliferation  of  nuclear  weapons  has  stimulated  interest  in 
seismic  discrimination  of  nuclear  weapons  tests  from  other  types  of  sources. 
Particularly  important  issues  for  CTBT  and  non-proliferation  monitoring  are 
extending  discrimination  capabilities  to  lower  magnitude  levels  and  into  new 
geographic  regions.  A  consequence  of  this  change  in  emphasis  is  that  different 
types  of  seismic  sources  begin  to  play  a  role  in  the  discrimination  problem.  One 
such  event  type  is  rockbursts  induced  by  mining  activity. 

Five  types  of  seismic  events  are  important  to  the  discrimination  problem: 
(1)  underground  nuclear  explosions,  (2)  earthquakes,  (3)  mine  or  quarry  blasts, 
(4)  construction  and  accidental  chemical  blasts,  and  (5)  rockbursts  or  tremors 
induced  by  mining.  Traditionally  seismic  discrimination  research  has  focussed 
on  the  first  two  of  these  source  types.  However,  as  the  threshold  of  interest  is 
pushed  to  lower  levels,  the  other  event  types  must  also  be  considered.  In  many 
areas  these  alternate  sources  are  frequent,  in  some  cases  they  may  be  large 
(particularly  for  rockbursts),  and  they  may  occur  anywhere  in  the  world. 
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Therefore,  a  reliable  scheme  for  monitoring  a  CTBT  worldwide  and  potential 
nuclear  weapons  testing  in  proliferant  countries  must  have  the  capability  to 
distinguish  these  different  event  types. 

1.2  Research  Objectives  and  Observations 

The  research  effort  described  here  has  been  directed  at  improving 
understanding  of  the  rockburst  or  mine  tremor  seismic  source  for  use  in 
discrimination.  Some  of  the  standard  teleseismic  discriminants  would  not  be 
expected  to  be  effective  for  rockbursts.  The  frequent  occurrence  of  these  kinds 
of  events  in  mining  areas  throughout  the  world  suggests  that  it  would  be  useful 
to  develop  techniques  to  facilitate  their  identification.  In  Figure  1  we  present  a 
map  of  the  locations  of  reported  rockbursts  or  mine  tremors,  which  were 
encountered  in  the  course  of  our  review  of  rockburst  activity  and  mechanisms.  It 
can  be  seen  that  these  events  have  been  reported  from  a  broad  range  of 
tectonic  regions  some  of  which  are  seismically  active  and  some  of  which  are 
not.  The  only  continent  which  appears  to  have  escaped  rockburst  activity  is 
South  America,  and  that  may  just  indicate  a  lack  of  reporting  or  failure  to 
investigate  this  kind  of  activity. 

For  purposes  of  this  study,  the  term  ■rockburst"  is  used  in  a  broad  sense 
to  include  any  type  of  stress-release  phenomenon  which  has  been  induced  by 
mining  activity  and  which  results  in  emission  of  seismic  signals.  Thus,  we  take 
rockbursts  to  include  mine  tremors,  collapses,  mine  bumps,  outbursts,  and  other 
similar  processes.  These  kinds  of  events  are  most  often  associated  with  mining 
in  rock  usually  at  substantial  depth;  but,  in  fact,  they  have  also  been  known  to 
occur  in  mining  near  or  even  at  the  earth's  surface.  During  the  mining  process, 
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Locations  of  mining  areas  with  reported  rockbursts  or  tremors 


potential  energy  is  added  to  the  surrounding  rock,  which  may  already  be 
prestressed  due  to  the  regional  tectonic  environment.  This  activity  causes  a 
readjustment  of  the  surrounding  stress  field  which  then  triggers  the  release  of 
seismic  energy.  Thus,  the  occurrence  of  rockbursts  in  mines  is  dependent  on 
the  mining  processes,  including  mine  volume,  rate  of  material  removal,  and 
support  conditions,  but  also  may  be  closely  related  to  the  regional  tectonic 
stress  and  strength  inhomogeneities  in  the  surrounding  rock. 

The  close  association  between  rockbursts  and  tectonic  stress  implies  that 
their  mechanisms  may  resemble  those  of  earthquakes.  Many  rockbursts  indeed 
do  appear  to  represent  release  of  tectonic  stress  on  preexisting  zones  of 
weakness  or  faults  near  the  mine  excavation  and,  therefore,  the  seismic  field 
which  they  generate  should  be  indistinct  from  that  of  an  earthquake  in  the  same 
region.  However,  some  rockbursts  are  more  complex  including  non-earthquake 
processes  such  as  rock  fall  or  pillar  failure.  Although  one  might  expect  that  the 
seismic  energy  release  in  these  latter  events  could  not  be  great,  we  found  that, 
in  fact,  one  of  the  largest  events  which  we  encountered  in  our  research  appears 
to  have  been  related  to  pillar  failure,  as  will  be  described  below. 

One  goal  of  mining  engineering  is  to  avoid  hazardous  conditions  in 
removing  ore,  so  we  generally  think  of  rockbursts  as  unusual  events.  However, 
in  nearly  all  mining  areas,  the  rockburst  hazard  has  not  been  completely 
eliminated;  and  in  some  areas  they  may  be  quite  frequent.  The  majority  of 
rockbursts  are  small  (mb  <  3);  however,  they  occasionally  (and  in  some  cases 
regularly)  exceed  5  mb-  In  comparison,  explosions  connected  with  quarrying, 
strip  mining,  or  construction  at  the  earth's  surface  seldom  have  magnitudes 
greater  than  3  -  4  mb;  and  explosions  used  in  development  of  underground 
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mines  are  usually  in  the  magnitude  range  1-3  mb-  Therefore,  in  many  parts  of 
the  world  mining-induced  events  may  pose  more  of  a  problem  to  discrimination 
monitoring  than  commercial  blasting  operations. 

The  ultimate  goal  of  discrimination  research  is  to  identify  techniques  to 
facilitate  distinction  of  underground  nuclear  explosion  tests  from  other  event 
types.  Some  features  of  rockbursts  suggest  that  they  may  behave  like 
earthquakes  in  many  aspects,  but  in  some  respects  the  behavior  should  be 
different.  Rockbursts  associated  with  induced  slip  along  faults  or  fractures 
should  behave  like  shallow  earthquakes,  and  earthquake  discrimination 
measures  should  be  useful  for  their  identification.  However,  rockbursts  are 
shallow  and  in  the  same  depth  range  as  nuclear  explosion  tests;  so  depth 
discriminants  will  not  generally  be  effective  in  separating  them.  Preliminary 
results  described  here  suggest  that,  with  respect  to  the  Ms  versus  mb 
discriminant  measure,  rockbursts  may  also  tend  to  look  explosion-like. 
However,  we  have  identified  features  in  the  spectra  of  teleseism ic  signals,  in  the 
relative  excitation  of  regional  phases,  and  in  the  spectra  of  regional  phases 
which  appear  to  offer  some  promise  for  distinguishing  rockbursts.  We  have  also 
found  evidence  of  complexity  in  the  first  motion  patterns  from  rockbursts  and  in 
the  signals  themselves  which  may  enable  their  distinction  from  nuclear 
explosion  tests.  Additional  observations  and  theoretical  studies  are  needed  to 
test  these  results. 

With  regard  to  controllability  of  rockbursts  and  the  potential  for 
clandestine  testing,  our  preliminary  finding  is  that  the  feasibility  of  predicting  the 
time  and  size  of  a  rockburst  event  in  a  mining  environment  could  provide  an 
opportunity  to  conceal  or  disguise  a  clandestine  nuclear  explosion  test.  There 
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is  evidence  that  rockbursts  can  be  predicted  and  deliberately  triggered  in  some 
mines.  Questions  regarding  the  maximum  size  of  events  which  may  be  induced 
and  the  precision  of  timing  predictions  are  still  outstanding  and  require  further 
study  to  reach  a  conclusion  on  this  scenario. 

1.3  Report  Organization 

The  report  is  organized  into  six  sections  including  this  introduction. 
Section  II  describes  the  seismic  databases,  including  rockbursts  and 
contrasting  source  types,  which  we  have  assembled  to  evaluate  rockburst 
characteristics  and  discrimination  issues.  The  discrimination  analyses  which 
have  been  performed  on  teleseismic,  far-regional  and  regional  waveforms  for 
Central  European  rockbursts  and  comparison  event  types  are  presented  in 
Section  III.  Section  IV  describes  similar  discriminant  analyses  on  the  South 
African  data.  In  Section  V  we  review  some  additional  discriminant  measures 
and  the  implications  of  rockburst  mechanisms  for  discrimination  including 
potential  problems  posed  by  rockburst  control  for  monitoring  clandestine 
testing.  Finally,  Section  VI  summarizes  our  observations  and  provides  some 
ideas  for  the  direction  of  future  investigations. 
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2.  Seismic  Database  for  Use  in  Rockburst 
Discrimination 

2.1  Seismic  Waveform  Database  Development 

In  our  previous  report  (Bennett  et  al.t  1993)  we  described  the 
characteristics  and  conditions  associated  with  rockbursts  in  different  areas 
around  the  world.  As  noted  there  and  in  the  introduction  to  this  report,  although 
some  rockbursts  or  associated  mine  tremors  are  large,  they  are  most  often  quite 
small  and  usually  not  very  well  recorded  at  normal  seismic  stations.  At  some 
mines  local  monitoring  networks  are  used  to  keep  track  of  mine  tremors  and 
warn  of  impending  hazards  to  mining  operations.  However,  in  many  countries 
there  appears  to  be  little  special  effort  to  account  for  rockburst  occurrences  and 
to  separate  them  from  other  types  of  natural  or  man-made  seismic  activity.  As  a 
result,  it  is  not  an  easy  matter  to  assemble  a  waveform  database  of  confirmed 
rockburst  events.  The  problem  is  further  complicated  by  the  fact  that  mines 
where  rockbursts  occur  may  also  use  chemical  explosions,  so  that  location  of 
an  event  at  a  mine  alone  does  not  distinguish  between  the  possible  source 
types. 

The  various  complexities  associated  with  database  compilation  led  us  to 
focus  on  two  areas  for  our  investigation  of  the  seismic  characteristics  of 
rockbursts:  Central  Europe  and  South  Africa.  In  both  of  these  source  regions, 
there  have  been  frequent  reports  of  rockbursts  or  mine  tremors.  In  each  area 
the  maximum  magnitudes  of  rockburst  events  have  exceeded  5.0  mb- 
Occurrence  of  large  events  permits  analysis  and  comparisons  of  teleseismic 
signals  from  each  source  region.  However,  in  each  area  there  also  exists  some 
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regional  monitoring  capability  which  should  enable  extending  the  source 
characterizations  to  smaller  events. 

For  both  Central  Europe  and  South  Africa,  it  was  not  our  goal  to  compile 
a  complete  database.  The  numerous  events  occurring  in  each  area  would 
make  such  a  task  impractical  for  a  research  project.  We  instead  sought  to  select 
representative  events  which  would  be  typical  of  rockbursts  from  each  region 
which  could  then  be  characterized  and  compared  with  other  types  of  sources. 
We  also  decided  to  consider  first  the  larger  events  with  reasonably  good  signal- 
to-noise  conditions  in  the  hope  that  characteristics,  which  might  stand  out  at  the 
higher  signal-to-noise  levels,  could  be  extended  to  the  smaller  events  where 
signals  might  be  weak. 

For  each  source  region  we  relied  on  digital  data  available  through  CSS, 
IRIS,  or  AFTAC.  Stations  from  which  rockburst  data  have  been  compiled 
include  SRO,  ASRO,  DWWSSN,  ARPA  regional  arrays,  Soviet  IRIS,  GSETT-2, 
and  AEDS  stations.  Events  which  are  included  in  the  compilation  produced 
good,  visible  signals  at  one  or  more  of  the  stations  whose  waveforms  we 
reviewed.  For  many  of  the  smaller  events,  signals  were  frequently  only  usable 
from  a  single  regional  station,  while  some  larger  rockbursts  produced  strong 
teleseismic  signals  at  ranges  of  100  degrees  or  more.  The  data  in  our  database 
usually  includes  all  the  available  short-period  and  broad-band  waveforms  with 
sampling  rates  of  20  per  second  or  higher.  For  regional  signals  we  tried  to 
retrieve  waveform  segments  corresponding  to  group  velocity  windows  from 
prior  to  P  through  Lg.  We  also  retrieved  long-period  surface-wave  data  at  one 
sample  per  second  for  many  events. 
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In  addition  to  the  rockbursts  from  each  source  region,  we  attempted  to 
include  in  our  database  other  source  types  which  could  be  used  for 
comparisons.  These  other  source  types  included  both  earthquakes  and 
chemical  blasts.  Wherever  possible  an  attempt  was  made  to  find  events  with 
similar  magnitudes  and  locations  to  the  rockbursts  for  use  in  the  comparisons. 
Because  coincident  sources  of  different  types  were  not  usually  available,  we 
alternatively  tried  to  pick  events  from  the  same  tectonic  region  and  at  about  the 
same  epicentral  distance  range.  We  also  compared  characteristics  of  seismic 
signals  between  similar  source  types  to  help  understand  variability  between 
events  of  a  common  type. 

As  we  described  in  our  previous  report  (cf.  Bennett  et  al.,  1993),  it  is  only 
rarely  that  published  information  exists  to  corroborate  events  as  rockbursts.  In 
most  cases  the  identification  of  events  as  rockbursts  is  based  on  inference 
derived  from  knowledge  that  the  events  occur  in  mining  areas  where  rockbursts 
have  been  reported  in  the  past  and  where  natural  seismicity  was  rare  or 
unknown  prior  to  the  initiation  of  mining  activity.  In  many  cases  the  mining 
areas  with  inferred  rockbursts  are  underground  so  that  blasting  is  either  not 
used  in  mine  development  or  it  involves  very  low  energy  release  (i.e.  smaller 
than  many  of  the  rockbursts).  As  a  result,  we  believe  it  is  unlikely  that 
mineblasts  or  natural  earthquakes  are  incorrectly  designated  as  rockbursts  for 
use  in  this  study.  Some  surface  quarrying  or  strip  mining  is  believed  to  use 
larger  blasts,  but  the  sources  of  such  events  are  usually  fairly  well-defined  and 
frequently  separate  from  the  rockburst  areas. 


9 


2.2  Central  European  Database 

Our  database  for  Central  Europe  consists  of  60  events  as  shown  in  Table 
1.  These  events  include  48  rockbursts  or  mine  tremors  (either  reported,  rb,  or 
probable  rockbursts,  prb),  7  earthquakes  (eqk),  4  mineblasts  (mb),  and  an 
accidental  underground  blast  (b).  A  few  small  events  had  no  reported 
magnitudes;  for  the  remainder  the  magnitudes  were  in  the  range  from  2.0  to  5.5 
mb-  The  majority  of  the  rockbursts  in  the  database  had  source  locations  in 
Poland  although  some  also  occurred  in  Germany.  The  earthquakes  had 
epicenters  in  the  Netherlands,  Germany,  Austria,  and  near  the  German/Czech 
border.  Mineblasts  had  locations  in  Austria  and  near  the  German/Czech  border 
and  the  Russian/Finnish  border.  The  accidental  underground  blast  occurred  in 
Switzerland.  The  map  in  Figure  2  shows  the  locations  of  these  60  events. 

The  Polish  rockbursts  are  clustered  in  two  areas  near  the  Czech  border: 
one  centered  near  51.1°N,  15.8°E,  and  the  other  near  51.0°N  19.0°E.  These 
correspond  to  the  copper  mining  area  near  Lubin  and  a  coal  mining  area  in 
Upper  Silesia  respectively.  Both  areas  have  had  long  histories  of  mining- 
induced  seismicity  (cf.  Gibowicz,  1984;  Bennett  et  al.,  1993).  The  largest 
German  rockburst  (5.4  mb)  has  attracted  considerable  attention  (cf.  Ahorner, 
1989;  Minkley,  1993)  and  appears  to  have  been  associated  with  the  extensive, 
nearly  simultaneous  collapse  of  multiple  pillars  in  a  mine  near  Volkershausen, 
as  will  be  discussed  more  completely  below.  The  large  earthquake  in  the 
Netherlands  (5.5  mb)  apparently  had  a  hypocenter  at  middle  to  lower  crustal 
depths  (see  below)  and  produced  numerous  aftershocks.  Little  is  known  about 
the  reported  mineblasts  in  the  database  except  that  they  were  generally  quite 
small  relative  to  the  other  events;  they  have  been  included  to  provide  another 
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Table  1.  Central 
Database 


Date 

Origin 

Lat  (N) 

Time 

08/1  5/80  '" 

20:03:04 

51.153 

08/25/80 

00:40:48 

51.115 

09/14/80 

15:39:13 

51.623 

09/30/80 

01:02:00 

50.350 

10/02/80 

17:23:44 

51 .319 

1  0/09/80 

04:46:59 

51.095 

11/29/80 

20:42:00 

151 .250 

12/21/80 

05:51:16 

i  51 .633 

12/23/80 

21:10:30 

151.636 

04/1 0/81 

04:59:47 

151.630 

07/1  2/81 

11:59:26 

1  50.522 

08/02/81 

03:25:16 

i  50.254 

08/1 1/81 

23:48:06 

149.861 

08/1  5/81 

05:40:14 

151.457 

08/1  9/81 

12:41:25 

152.103 

09/23/81 

23:53:41 

151.153 

11/03/81 

05:42:44 

151.027 

11/2  5/81 

10:14:20 

151.055 

12/1  1/81 

02:02:32 

!  51 .161 

01/1  6/82 

11:15:59 

[51.610 

01/22/82 

03:09:03 

151.156 

03/19/82 

14:48:05 

1  50.988 

04/12/82 

12:45:37 

04/1  8/82 

51 .651 

04/22/82 

21:22:38 

51.153 

05/23/82 

18:22:41 

50.851 

06/04/82 

10:44:34 

50.537 

06/08/82 

51.444 

06/08/82 

16:33:40 

51.156 

06/1  3/82 

21:15:27 

51 .689 

06/15/82 

20:00:59 

51.620 

mamm 

09:30:42 

49.915 

06/1  9/82 

23:36:29 

QEiEm 

07/1  4/82 

00:56:37 

51.681 

08/24/85 

06:08:40 

50.275 

10/25/86 

20:07:23 

51 .640 

05/31/87 

02:54:27 

51 .600 

06/20/87 

51.480 

03/1  3/89 

13:12:15 

50.710 

03/2  5/89 

12:46:40 

51.620 

Events  In  Current 


Lon  (E) 

Mag 

ID 

15.892 

3.6 

orb 

15.745 

ITU  H  I— 

16.200 

3.6 

prb 

18.910 

- 

prb 

15.571 

3.9 

prb 

15.687 

prb 

19.400 

■ 

prb 

16.259 

4.2 

prb 

16.375 

3.7 

prb 

16.180 

4.3 

prb 

19.012 

4.1 

prb 

18.779 

- 

prb 

18.425 

4.1 

prb 

15.954 

4.0 

prb 

17.574 

3.4 

Prb . 

15.818 

2.8 

JL£b _ 

15.867 

twwmmm 

prb 

15.828 

3.9 

prb 

15.817 

3.9 

Prb . 

16.281 

4.2 

prb 

15.516 

3.4 

Prb  . 

15.155 

3.3 

prb 

1 

15.966 

prb 

1 

16.085 

3.5 

prb 

15.860 

Elg— m 

prb 

15.158 

19.062 

DH 

prb 

15.961 

4.0 

prb 

15.831 

3.4 

prb 

15.946 

4.1 

prb 

16.278 

3.7 

Elb 

18.470 

3.8 

prb 

1 

15.986 

3.7 

prb 

16.098 

4.0 

prb 

7.917 

4.3 

eqk 

16.080 

4.9 

prb 

16.270 

4.9 

prb 

16.050 

4.7 

prb 

9.900 

rb 

16.090 

prb 

11 


Date 


Origin 

Time 


Lat  (N) 

Lon  (E) 

51 .438 

16.120 

51  .600 

16.010 

60.914 

29.131 

51 .542 

16.227 

47.962 

16.204 

52.281 

7.761 

49.288 

6.930 

50.312 

12.441 

50.260 

19.220 

51  .583 

16.085 

50.140 

12.850 

50.183 

12.715 

51 .528 

16.392 

51.561 

16.115 

51.374 

6.133 

51 .043 

5.910 

51 .099 

5.885 

50.748 

6.039 

47.500 

14.900 

46.730 

8.330 

06/29/89 


10/28/89 


1  1/29/90 


04/30/91 


05/02/91 


05/1  6/91 


05/1  6/91 


05/1  9/91 


05/21/91 


lEHEFfEU 


05/26/91 


05/26/91 


05/28/91 


1  0/22/91 


04/13/92 


04/1  3/92 


04/1  3/92 


04/14/92 


07/03/92 


1  1/02/92 


17:33:51 


01:35:53 


12:29:1 1 


03:40:39 


10:15:20 


02:06:17 


10:44:58 


03:22:13 


16:49:09 


1 1:01:01 


19:42:55 


03:52:50 


19:19:25 


01:20:05 


03:49:41 


06:02:10 


01:06:46 


07:53:01 


15:13:26 


113 


67  N 
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event  locations  for  the  Central  European  database 


source  type  for  signal  comparisons.  The  accidental  explosion  in  Switzerland 
had  a  magnitude  of  4.2  mb  and  was  the  result  of  the  detonation  of  ammunition 
in  a  storage  bunker.  This  event  appears  to  be  the  largest  blast  in  the  recent 
history  of  Central  Europe  which  has  been  well-recorded  on  high-quality  seismic 
stations  and  thus  provides  another  valuable  comparison;  however,  little 
information  about  the  source  of  this  event  has  been  published  at  this  time. 

In  a  previous  report,  Bennett  et  al.  (1993)  reviewed  the  general 
characteristics  of  seismic  signals  from  events  in  Central  Europe.  The  largest 
events  produce  strong  body  wave  signals  at  teleseismic  stations.  However, 
many  of  the  small  rockbursts  and  mineblasts  are  not  well  recorded  at 
teleseismic  stations;  and,  therefore,  regional  signals  are  important  to  their 
identification.  In  this  study  we  have  performed  discrimination  analyses  on  data 
for  several  large  Central  European  events  recorded  at  teleseismic  stations, 
including  AEDS  network  stations.  We  have  also  done  comparative  analyses  of 
different  Central  European  event  types  at  nearly  common  far-regional  epicentral 
distances  recorded  at  the  ARPA  regional  array  at  NORESS.  Finally,  we  have 
conducted  more  detailed  discrimination  analyses  of  several  different  types  of 
smaller  events  recorded  at  about  the  same  smaller  regional  distances,  but 
different  azimuths,  from  the  GDSN  station  GRFO. 

2.3  South  African  Database 

Our  database  for  South  Africa  includes  79  events  (cf.  Table  2).  Except  for 
a  few  events  of  unknown  magnitude,  the  magnitudes  are  in  the  range  from  2.3 
to  5.5  Ml.  The  majority  of  the  events  (61)  are  thought  to  be  rockbursts  or  mine 
tremors  (prb)  based  on  their  location  within  the  area  of  deep  gold  mining  in  the 
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Table  2.  South  African  Events  in  Current  Database 


Date 


Origin 

Time 


Lat  (S)  |  Lon  (E) 


01/28/80 


02/09/80 


02/17/80 


03/22/80 


04/03/80 


04/08/80 


05/06/80 


06/12/80 


06/13/80 


01/28/81 


02/18/81 


07/07/81 


11/05/81 


1  1/19/81 


1  1/22/81 


1  1/22/81 


11/26/81 


12/01/81 


12/01/81 


12/14/81 


12/15/81 


12/15/81 


03/26/82 


03/28/82 


04/01/82 


04/02/82 


04/02/82 


04/09/82 


04/13/82 


04/13/82 


05/09/82 


06/07/82 


06/16/82 


HI 


06/27/82 


06/28/82 


09/01/82 


SEBEEai 


11/12/82 


12/02/82 


06:30:57 


13:51:10 


21:01:22 


02:14:10 


22:45:38 


00:39:35 


21:07:25 


03:03:42 


21:15:03 


16:20:33 


08:28:20 


12:52:42 


20:19:31 


17:19:42 


03:31:27 


05:15:59 


07:16:49 


09:04:34 


15:11:45 


14:41:01 


14:11:15 


13:41:24 


15:51:37 


07:11:01 


12:57:09 


20:08:54 


01:52:37 


1 1:26:00 


11:26:51 


07:07:02 


■ram 


17:57:22 


04:18:52 


00:36:36 


09:20:08 


11:38:20 


09:24:21 


06:11:32 


20:45:34 


26.396 


27.802 


27.627 


26.120 


26.081 


26.293 


26.926 


26.991 


26.845 


26.073 


26.625 


26.552 


29.950 


28.000 


26.400 


26.380 


24.080 


25.207 


26.400 


26.320 


27.600 


26.360 


27.660 


26.270 


30.040 


26.830 


26.220 


26.800 


27.920 


27.900 


29.600 


27.900 


23.500 


26.140 


26.760 


26.880 


27.930 


26.200 


26.906 


30.600 


27.463 


26.235 


26.865 


27.752 


27.703 


27.653 


26.998 


26.991 


26.769 


27.611 


26.607 


27.026 


27.370 


26.840 


27.500 


27.570 


30.980 


27.912 


27.370 


27.330 


27.100 


27.530 


31.100 


28.220 


19.400 


26.750 


28.090 


26.600 


26.780 


26.800 


27.060 


26.800 


26.100 


27.710 


26.540 


26.810 


26.830 


29.900 


26.750 
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Date 


Origin 

Time 

Lat  (S) 

Lon  (E) 

Mag 

ID 

22:03:04 

26.900 

26.600 

22:03:59 

26.830 

26.720 

16:26:41 

imUTT'l 

28.493 

1 1:49:51 

HOI 

26.720 

13:09:57 

27.990 

26.800 

00:55:22 

26.200 

28.130 

16:13:59 

26.220 

28.150 

10:48:51 

26.890 

26.660 

13:38:07 

26.140 

27.890 

31.190 

24.250 

16:58:21 

26.220 

28.210 

06:05:48 

33.390 

19.270 

21:19:26 

30.830 

21  .500 

26.350 

27.320 

00:33:36 

29.470 

25.020 

00:33:43 

29.200 

24.800 

03:05:34 

29.540 

24.900 

08:59:12 

28.020 

26.880 

04:43:13 

E-KEHMK 

27.700 

23:16:47 

30.060 

25.790 

07:38:24 

26.421 

27.387 

14:40:02 

BEH 

26.650 

16:47:58 

26.190 

27.850 

06:45:07 

20.140 

26.250 

21:23:10 

26.800 

26.520 

15:31:52 

29.580 

26.710 

16:00:50 

26.780 

26.600 

20:45:02 

27.940 

26.750 

04:59:10 

26.920 

26.570 

07:06:31 

26.270 

27.430 

15:04:21 

26.980 

26.680 

20:41:59 

26.304 

27.486 

06:41:17 

26.807 

26.639 

10:14:33 

27.985 

26.734 

16:17:10 

26.872 

26.736 

23:08:24 

28.014 

26.727 

00:31:39 

28.122 

26.789 

18:28:55 

26.895 

IIITT'BU 

12/11/82 


02/22/83 


05/24/83 


05/24/83 


06/01/83 


06/03/83 


06/06/83 


07/17/83 


07/31/83 


08/01/83 


08/04/83 


08/26/83 


08/26/83 


09/05/83 


09/05/83 


09/09/83 


09/29/83 


10/01/83 


1  1/02/83 


12/27/83 


01/28/84 


03/21/84 


06/04/84 


08/1 1/84 


08/19/84 


01/01/86 


02/10/86 


08/1  1/86 


09/15/86 


10/28/86 


12/29/86 


01/05/88 


01/25/89 


02/27/90 


09/26/90 


1  1/14/90 


11/03/91 


12/01/91 


26.664 


Witwatersrand  basin  near  27°S  27°E.  The  remaining  18  events  have  locations 
outside  this  area  and  are,  therefore,  classified  here  as  probable  earthquakes 
(eqk),  although  in  only  a  few  cases  does  actual  corroboration  of  the  source  type 
exist  in  the  ISC  Bulletin  or  South  African  Seismicity  Bulletins.  The  map  in 
Figure  3  shows  the  locations  of  the  events  in  our  South  African  database. 

In  our  previous  report  (Bennett  et  al.,  1993)  we  described  some  of  the 
general  characteristics  of  rockburst  occurrence  in  the  South  African  gold  mining 
district.  The  induced  activity  there  follows  a  somewhat  arcuate  trend  coincident 
with  the  mine  locations.  There  has  been  an  abundance  of  reports  in  the 
scientific  literature  on  the  occurrence  and  mechanisms  of  South  African 
rockbursts.  The  events  are  generally  deep  and  frequent,  with  nearly  85  events 
per  year  in  the  magnitude  range  3.4  to  5.9  Ml  during  the  period  1971  to  1979. 
Because  of  the  frequency  and  long  history  of  induced  seismicity  and  extensive 
monitoring  program  in  South  Africa,  more  is  probably  known  about  the 
mechanisms  and  relation  of  their  occurrence  to  excavation  operations  there 
than  in  any  other  area.  We  will  summarize  the  general  state  of  understanding  of 
South  African  rockburst  mechanisms  in  a  subsequent  section  of  this  report. 

Bennett  et  al.  (1993)  noted  that  the  best  available  GDSN  data  for  South 
African  rockbursts  came  from  the  regional  station  at  SLR  and  from  a  few 
teleseismic  stations,  including  BCAO  in  Central  Africa  and  ZOBO  in  Bolivia, 
South  America.  In  this  report  we  have  added  valuable  data  from  teleseismic 
AEDS  network  stations  for  several  South  African  events.  We  will  describe  some 
parametric  observations  as  well  as  spectral  characteristics  of  the  teleseismic 
AEDS  station  signals  for  these  events.  We  have  also  performed  additional 
discrimination  analyses  of  the  BCAO  long-period  data  for  several  of  the  larger 
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Figure  3.  Map  of  event  locations  for  the  South  African  database. 


events.  As  was  the  case  for  the  larger  Central  European  events  at  GRFO,  large 
South  African  mine  tremors  frequently  exceed  the  recording  threshold  at  the 
nearby  SLR  station.  Therefore,  the  near-regional  monitoring  effort  for  South 
African  rockbursts  must  focus  on  smaller  events.  In  a  latter  section  we  will 
summarize  some  of  the  characteristics  of  near  regional  signals  from  South 
African  rockbursts  based  on  the  SLR  observations  and  their  relevance  to 
discrimination. 
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3.  Central  European  Events 

3.1  Teleseismic  Signals 

We  begin  our  analyses  of  Central  European  rockbursts  with  a 
comparison  of  teleseismic  signals  for  three  large  Central  European  events 
recorded  at  AEDS  stations.  The  events  used  in  this  comparison  were  the 
4/13/92  Netherlands  earthquake  {5.5  mb),  the  3/13/89  German  rockburst  (5.4 
mb),  and  the  6/29/89  Polish  rockburst  (4.7  mb).  After  reviewing  the  available 
waveform  data,  we  selected  the  records  from  eight  common  AEDS  stations  for 
each  event.  The  stations  were  at  distances  in  the  range  50°  <  A  <  80°. 

With  regard  to  general  teleseismic  waveform  features,  the  initial  P 
waveforms  for  the  rockbursts  at  common  stations  appear  somewhat  more 
complex  than  for  the  first  earthquake  P  signals.  At  least  two  factors  may 
contribute  to  the  observed  differences.  Since  the  rockbursts  are  shallower  than 
the  earthquakes,  site  response  of  the  shallow  near-source  crust  may  produce 
longer  duration  and  complexity  compared  to  the  earthquake.  A  second  factor  is 
complexity  of  the  mechanism.  Some  rockburst  sources  appear  to  be  more 
complex  including  multiple  events  or  episodes.  In  fact,  comparison  of  the 
3/13/89  rockburst  P-wave  signals  with  the  6/29/89  rockburst  signals  suggests 
that  some  rockbursts  are  more  complex  than  others.  The  3/13/89  rockburst 
appears  to  have  a  precursory  phase  arriving  about  two  seconds  prior  to  the 
main  P  at  all  the  teleseismic  AEDS  stations.  This  phase  is  approximately  one- 
third  the  amplitude  of  the  main  P  phase.  There  are  also  indications  in  the 
rockburst  record  of  two  or  more  phases  occurring  after  the  main  P  within  5  -  6 
seconds.  Because  the  intervals  of  these  multiple  P  phases  appear  to  be  nearly 
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the  same  at  all  stations  (i.e.  no  evidence  of  moveout),  it  seems  unlikely  that  they 
are  caused  by  multiple  mantle  phases  but  are  more  probably  produced  by 
multiple  events  at  the  source.  In  fact,  such  complexity  in  the  source  of  the 
3/13/89  rockburst  may  explain  some  questions  involved  in  interpreting  the 
mechanism  from  P-wave  first  motions,  as  will  be  discussed  in  a  subsequent 
section. 

One  other  difference  in  the  earthquake  signals  is  the  appearance  of  a 
rather  strong  depth  phase.  This  phase  arrives  about  7.5  -  8  seconds  after  the 
initial  P  and  is  quite  strong  on  the  records  at  nearly  all  the  AEDS  stations. 
Assuming  this  phase  to  be  pP  the  difference  in  time  relative  to  P  would  indicate 
a  source  depth  of  about  20  km.  This  agrees  with  the  source  depth  of  21  km 
reported  by  NEIC  for  the  4/13/92  Netherlands  earthquake. 

3.2  Spectral  Comparisons  of  Teleseismic  P 

To  develop  a  better  understanding  of  possible  spectral  differences  and 
variability  in  the  P  wave  signals  from  rockbursts  and  other  sources,  we 
performed  spectral  analyses  on  the  P  signals  at  the  AEDS  stations  for  the 
4/13/92  earthquake,  the  3/13/89  rockburst,  and  the  6/29/89  rockburst.  We 
analyze  the  results  here  in  terms  of  the  P/P  spectral  ratios  between  events  as 
recorded  at  common  stations.  Thus,  we  represent  the  observed  P-wave 
spectrum  at  a  station  as 

A(o>)  =  S0{co)  Pr(co)  Ts(co)  (1 ) 

where  So(a))  is  the  source  spectrum  (including  source  mechanism  and  source 
site  response  effects),  Pr(cd)  is  a  propagation  path  response,  and  Ts(co)  is  the 
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station  site  and  instrument  transfer  function.  Then,  for  two  different  events  the 
P/P  spectral  ratio  is: 


R  a  A,(a>)  =  SJco)  PB1(m)  Tsi(m) 

%  A»  So2(co)  PR2(oj)  TS2((b)  (2> 

Since  we  are  taking  the  ratios  at  common  stations,  Tsi  =  TS2  and  we  will  also 
assume  that  Pri  «  PR2  since  all  sources  are  in  the  same  region  and,  therefore, 
have  about  the  same  transmission  path  to  the  common  stations.  Thus,  the 
spectral  ratio, 


v*>  *  ftS 


(3) 


should  provide  an  indication  of  source  spectral  differences  between  events. 

For  the  Central  European  earthquake  and  rockbursts,  we  selected  P- 
wave  window  segments  7.5  seconds  long  from  each  of  eight  teleseismic, 
common  AEDS  stations  and  a  noise  segment  of  similar  duration  from  the  time 
segment  preceding  P.  Fourier  spectra  were  computed  for  both  the  P  and  the 
noise  segments  corresponding  to  each  waveform.  Spectral  bands  in  which 
signal-to-noise  (S/N)  ratios  were  strongest  for  each  station  and  event  were 
identified.  For  each  common  station  we  then  computed  the  P/P  spectral  ratios 
for  the  3/13/89  rockburst  to  the  4/13/92  earthquake  and  for  the  6/29/89  rockburst 
to  the  4/13/92  earthquake.  The  ratios  were  computed  over  the  frequency  band 
in  which  the  S/N  ratio  was  good  {generally  greater  than  one)  at  each  station. 
The  average  P/P  ratios  over  all  the  common  stations  were  next  computed  using 
values  of  the  ratios  for  which  S/N  was  good.  This  last  averaging  process  should 
further  reduce  any  station  effects  which  might  persist  in  the  ratios. 

The  average  ratios  with  plus-  and  minus-one-sigma  (±  1  o)  bounds  about 
the  mean  ratios  can  then  be  compared.  These  are  presented  in  Figure  4.  At 
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Figure  4.  Average  P/P  spectral  ratios  and  ±1a  bounds  determined  from  teleseismic  AEDS 
stations  for  large  Central  European  rockbursts  and  the  large  Netherlands  earthquake. 


low  frequencies  (0.5  -  1.5  Hz)  the  P/P  spectral  ratios  for  3/13/89 
rockburst/4/13/92  earthquake  is  relatively  flat  and  oscillates  slightly  around  a 
value  near  1.0.  This  is  consistent  with  the  comparable  mb  values  of  the  two 
events.  The  average  ratio  has  a  prominent  peak  near  2  Hz  indicating 
enhancement  of  the  rockburst  P  signal  relative  to  the  earthquake,  or  deficiency 
of  the  earthquake  relative  to  the  rockburst,  within  that  band.  Ignoring  the  peak 
near  2  Hz,  the  general  behavior  of  the  P/P  spectral  ratio  appears  to  represent  a 
rather  steady  decline  toward  higher  frequencies,  suggesting  that  the 
earthquake  has  more  high-frequency  P-wave  energy  than  the  rockburst. 

For  the  6/29/89  rockburst  comparison  to  the  4/13/92  earthquake  shown  at 
the  right  of  Figure  4,  the  P/P  spectral  ratios  are  generally  lower.  This  is  probably 
primarily  indicative  of  the  smaller  magnitude  of  this  rockburst  event  relative  to 
the  earthquake.  The  spectral  ratio  again  shows  some  rather  prominent  peaks 
near  2  Hz  and  also  near  0.8  Hz.  Unlike  the  other  example,  however,  the  trend 
over  the  broader  frequency  band  is  generally  upward  toward  higher 
frequencies.  Several  factors  could  be  responsible  for  this  difference  in  the 
spectral  ratios.  The  increase  with  frequency  observed  for  the  6/29/89  event 
might  be  explained  partly  by  a  shift  toward  higher  corner  frequency  for  the 
smaller  event.  An  alternative  explanation,  which  is  probably  most  interesting 
from  the  standpoint  of  discrimination,  is  that  the  difference  in  mechanism 
between  the  German  event  (thought  to  have  been  a  collapse)  and  the  Polish 
rockburst  (thought  to  have  been  more  tectonic  in  character)  is  the  source  of  the 
spectral  differences.  Two  more-mundane  explanations  for  the  observed 
differences  in  behavior  might  be  differences  in  the  near-source  site  response 
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between  the  two  rockbursts  or  some  degree  of  noise  contamination  at  higher 
frequencies  in  the  spectrum  for  the  smaller  event. 

Thus,  the  results  of  the  teleseismic  P-wave  analyses  for  the  Central 
European  rockbursts  are  not  so  clear.  There  may  be  variability  between  events 
with  some  types  of  rockbursts  (tectonic)  tending  to  be  enriched  in  higher 
frequencies  compared  to  collapse  type  events.  However,  additional  work  is 
needed  to  resolve  alternative  explanations  and  tie-down  source  identification 
procedures.  This  work  should  include  analyses  of  additional  observations  as 
well  as  theoretical  investigations  of  variability  in  source  mechanisms. 

3.3  Far-Regional  Signal  Analyses 

In  addition  to  the  teleseismic  data,  regional  signals  are  expected  to 
provide  a  valuable  contribution  to  rockburst  identification,  particularly  as  events 
approach  smaller  magnitude  levels.  As  part  of  our  investigation  of  rockburst 
signal  characteristics,  we  analyzed  the  behavior  of  several  events,  from  the 
database  described  in  Section  II,  recorded  at  far-regional  stations.  These 
analyses  again  focus  on  different  event  types  recorded  at  common  stations. 

Figure  5  shows  a  map  with  the  locations  of  three  fairly  large  rockbursts 
and  the  large  Netherlands  earthquake  in  Central  Europe  along  with  the 
locations  of  three  good  far-regional  stations  which  would  be  expected  to 
provide  good  signals  from  the  events.  Figure  6  shows  a  comparison  of  the 
vertical-component  signals  recorded  at  these  stations  from  the  3/13/89  German 
rockburst  (5.4  mb)  and  the  4/13/92  Netherlands  earthquake  (5.5  mb).  For  the 
two  events  the  distances  to  the  individual  stations  are  roughly  comparable. 
Several  differences  are  apparent  in  the  time-domain  signal  behavior.  At  NRAO 
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Figure  5.  Locations  of  several  rockbursts,  an  earthquake,  and  a  mineblast  in  Central  Europe 
used  in  discrimination  analysis.  Inverted  triangles  show  selected  far-regional  stations. 
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Figure  6.  Comparison  of  vertical-component  waveforms  at  three  regional  stations  from  the  1989  German 
rockburst  and  the  1992  German/Netherlands  border  earthquake. 


the  Lg  and  P  signals  appear  enhanced  relative  to  S  for  the  rockburst  on  these 
broader-band  recordings.  At  OBN  and  TOL  S  appears  somewhat  more 
prominent  for  the  earthquake  than  for  the  rockburst,  but  there  is  also  a  very 
sharp  P  phase  in  the  earthquake  record  at  OBN.  At  both  OBN  and  TOL  there  is 
evidence  of  significant  long-period  surface-wave  excitation  for  both  events. 

We  performed  a  band-pass  filter  analysis  of  the  NORESS  records  for  the 
earthquake  and  the  three  rockbursts  from  Figure  5.  Each  record  was  passed 
through  a  set  of  overlapping  band-pass  filters  with  similar  characteristics  but 
different  cut-off  frequencies.  Comparing  the  1989  German  rockburst  and  the 
1992  Netherlands  earthquake  in  Figure  7,  there  is  an  apparent  enhancement  of 
P  and  Lg  relative  to  S  in  the  rockburst  compared  to  the  earthquake  which 
persists  at  all  frequencies.  There  also  appears  to  be  a  tendency  in  the  rockburst 
P  to  be  more  prominent  relative  to  other  phases  at  high  frequencies.  This 
appears  to  be  contrary  to  our  teleseismic  observations  which  indicated  more 
high-frequency  energy  in  the  earthquake  P  than  in  the  rockburst  P  for  these 
same  events,  but  the  prominence  of  P  here  may  be  related  more  to  a  more  rapid 
decline  in  S  and  other  phases  for  the  rockburst  than  to  actual  P  enhancement. 
Spectral  analyses  of  the  regional  P  signals  for  these  events  is  discussed  in 
greater  detail  later  in  this  section.  The  same  band-pass  filter  analyses  of  the 
other  two  rockbursts  are  shown  in  Figure  8.  In  these  examples  we  do  not  see 
the  strong  Lg  relative  to  S  at  all  frequencies;  only  in  the  1989  Polish  rockburst  is 
the  Lg  prominent  and  primarily  at  lower  frequencies.  However,  the  P  signals  in 
these  rockbursts  again  appear  to  be  stronger  relative  to  S  than  they  were  for  the 
earthquake  in  Figure  7.  Once  again  there  is  a  tendency  for  the  rockburst  P  to 
increase  in  prominence  at  high  frequencies  relative  to  other  phases. 
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989  German  rockburst  and  1992  German/Netherlands  border  earthquake. 
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German  rockburst  and  1989  Polish  rockburst. 


Figure  9  presents  the  results  from  Figures  7  and  8  in  terms  of  parametric 
measurements  of  the  peak  S/P  amplitude  ratios  in  each  frequency  band.  The 
S/P  ratios  for  the  1992  Netherlands  earthquake  are  in  the  3  -  5  range  declining 
only  slightly  with  frequency.  S/P  ratios  for  the  rockbursts  are  generally  in  the  1  - 
3  range  and  show  more  decline  with  frequency.  Similar  analyses  including  a 
few  additional  events  from  the  database  continued  to  suggest  a  tendency  for  the 
rockbursts  and  earthquakes  from  Central  Europe  recorded  at  NORESS  to 
separate  on  the  basis  of  S/P  ratio.  However,  it  should  be  noted  that  all  events 
used  in  these  analyses  do  not  have  the  same  propagation  path  to  NORESS; 
and  we  have  not  attempted  here  to  adjust  for  propagation  differences.  This  is 
left  to  future  refinement  of  the  method. 

3.4  Further  Investigation  of  Far-Regional  P  Spectra 

To  better  understand  the  behavior  of  the  far-regional  P-wave  spectra 
observed  at  NORESS,  we  analyzed  the  signals  from  the  3/13/89  German 
rockburst  and  the  4/13/92  Netherlands  earthquake  already  described  plus  the 
10/28/89  Polish  rockburst  (4.7  mb)  and  a  mineblast  near  the  Russian/Finnish 
border  (1 1/29/90,  2.7  Ml)  as  shown  in  Figure  5  above.  All  four  events  are 
located  at  comparable  epicentral  distances  from  NORESS  - 1120  km,  1135  km, 
1055  km,  and  950  km  respectively.  All  four  events  produced  good  P  signals  at 
NORESS.  Signal  levels  were  well  above  noise  out  to  frequencies  of  10  Hz  or 
more,  although  the  P  signal  for  the  smaller  magnitude  mineblast  appears  to  be 
contaminated  by  noise  at  low  frequencies  (f  <  2  Hz). 

We  analyzed  the  spectral  behavior  in  terms  of  P/P  spectral  ratios  which 
are  shown  in  Figure  10.  For  the  3/13/89  rockburst/4/13/92  earthquake 
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Filter  Passband 

Figure  9.  Comparison  of  S/P  ratios  at  NORESS  for  various  filter  passbands  for  selected 
Central  European  rockbursts  and  Netherlands/German  border  earthquake. 
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Figure  10.  P/P  spectral  ratios  determined  from  far-regional  signals  at  NORESS  for  large 
Central  European  rockbursts  and  the  large  Netherlands  earthquake. 


comparison,  the  P/P  ratio  is  near  a  value  of  one  in  the  vicinity  of  1  Hz.  This  is 
consistent  with  the  similar  magnitude  levels  of  the  two  events  and  with  our 
teleseismic  results  presented  earlier  in  Section  3.2.  The  ratio  has  a  peak  near  2 
Hz,  again  consistent  with  the  teleseismic  results,  but  generally  falls  off  rather 
steadily  to  a  value  less  than  0.1  at  10  Hz.  The  P/P  ratio  for  the  10/28/89 
rockburst/4/13/92  earthquake  is  somewhat  less  well-behaved;  the  ratio  shows 
some  spectral  holes  which  produce  rapid  fluctuations  in  level.  The  ratio  is  well 
below  1.0  near  1  Hz  indicative  of  the  smaller  magnitude  of  the  Polish  rockburst 
relative  to  the  earthquake.  The  P/P  spectral  ratio  tends  to  decrease  with 
frequency  out  to  about  5  Hz  and  then  levels  off  up  to  10  Hz.  Reasons  for  the 
leveling  off  of  the  ratio  are  not  clear  but  could  be  caused  by  noise  contamination 
or  a  shift  in  corner  frequency  for  the  smaller  event  or  also  possibly  by 
attenuation  differences  between  the  two  propagation  paths. 

As  another  check  on  the  far-regional  P  signal  behavior,  we  looked  at  the 
signals  from  the  3/13/89  German  rockburst  and  the  4/13/92  Netherlands 
earthquake  recorded  at  the  AEDS  station  in  Sonseca,  Spain.  The  epicentral 
distances  in  this  case  were  1310  km  and  1285  km  from  the  German  rockburst 
and  the  earthquake  respectively.  The  S/N  comparisons  for  the  earthquake 
indicate  good,  strong  signals  which  are  well  above  noise  over  a  broad  band 
from  low  frequency  out  to  beyond  5  Hz.  However,  for  the  rockburst  the  signal  is 
unusually  low,  approaching  noise  near  3  -  5  Hz.  We,  nevertheless,  computed 
the  P/P  spectral  ratio  over  the  band  from  0.5  to  5  Hz  as  shown  in  Figure  1 1 .  The 
ratio  is  generally  very  low  (less  than  0.1)  over  most  of  the  band  and  shows  a 
tendency  to  decline  with  increasing  frequency.  We  assume  there  is  some 
problem  in  the  calibration  information  provided  for  the  Sonseca  station  for  the 
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Figure  1 1 .  P/P  spectral  ratio  determined  from  far-regional  signals  at  the 
Sonseca,  Spain  AEDS  station  for  the  large  Central  European  rockburst 
and  earthquake. 


3/13/89  rockburst  waveforms.  This  assumption  appears  to  be  supported  by 

observations  of  these  same  two  events  on  the  nearby  GDSN  station  at  TOL, 

where  the  amplitude  difference  between  the  two  events  is  seen  to  be  only  about 
30  percent.  Disregarding  the  difference  in  level,  the  trend  in  the  ratio  in  Figure 
1 1  again  would  indicate  a  more  rapid  decline  in  high-frequency  energy  in  the 
3/13/89  rockburst  than  in  the  4/13/92  earthquake. 

For  the  11/29/90  mineblast/4/13/92  earthquake  comparison  the  P/P 
spectral  ratio  is  presented  in  Figure  12.  We  see  that  the  ratio  has  a  low  value 
much  less  than  0.1  -  near  1  -  2  Hz  in  keeping  with  the  much  lower  magnitude  of 
the  mineblast  compared  to  the  earthquake.  However,  the  ratio  shows  a 
tendency  to  increase  with  frequency  out  to  10  Hz.  The  behavior  might  represent 
some  combination  of  the  excitation  mechanism  and  corner  frequency.  Scaling 
analyses  would  be  useful  to  identify  to  what  extent  the  higher  corner  frequency 
of  the  smaller  magnitude  event  (viz.  mineblast)  influences  the  spectral  ratio. 
Higher  frequencies  in  the  mineblast  might  also  be  attributed  to  propagation  path 
differences  if  it  can  be  assumed  that  the  path  from  NORESS  into  Finland  has 
lower  P-wave  attenuation  than  the  path  to  the  south  into  Europe  and  the 
Netherlands. 

We  would  conclude  from  these  observations  that  the  far-regional  P-wave 
observations  for  Central  Europe  are  generally  consistent  with  the  teleseism ic  P- 
wave  observations  reported  above.  Both  teleseismic  and  far-regional  P  waves 
for  the  3/13/89  rockburst  appear  to  have  less  high  frequency  energy  than  similar 
P  waves  for  the  4/13/92  earthquake.  However,  interpretation  of  the  behavior  for 
the  other  events  is  more  problematic.  The  Polish  rockburst  (10/28/89),  for  which 
we  analyzed  the  far-regional  P  signals,  appears  to  have  less  high-frequency 
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Figure  1 2.  P/P  spectral  ratio  determined  from  far-regional  signals  at 
NORESS  for  a  mineblast  near  the  Finnish/Russian  border  and  a  large 
Central  European  earthquake. 
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energy  than  the  earthquake;  while  the  Polish  rockburst  (6/29/89),  for  which  we 
analyzed  teleseismic  P,  shows  the  opposite  trend.  A  reasonable  interpretation 
of  the  results  appears  to  be  that  rockburst  or  mine  tremor  mechanisms  may  vary 
in  the  efficiency  with  which  they  excite  P-wave  energy  in  different  frequency 
bands.  Theoretical  model  studies  and  improved  knowledge  of  path  attenuation 
effects  could  provide  important  evidence  to  test  this  hypothesis. 

3.5  Nearer-Regional  Observations  - 

Our  final  analyses  of  the  Central  European  events  focuses  on  nearer- 
regional  observations  at  station  GRFO.  In  these  analyses  we  primarily  relied  on 
data  from  three  events:  the  8/15/81  Polish  rockburst  (4.0  Ml),  the  4/14/92 
aftershock  (4.3  Ml)  of  the  Netherlands  earthquake,  and  the  11/02/92  Swiss 
munitions  blast  (4.0  Ml).  As  noted,  the  events  all  had  roughly  the  same 
magnitudes;  and,  as  can  be  seen  in  the  map  in  Figure  13,  they  are  each  located 
at  about  the  same  epicentral  distance  from  GRFO  (390  km,  387  km,  and  393  km 
respectively). 

Figure  14  shows  the  vertical  component  records  for  the  three  events. 
Approximate  times  of  Pn,  Pg,  and  Lg  phases  are  identified  on  the  records.  The 
Lg  phases  are  the  largest  on  the  broadband  recordings  for  all  three  events.  S 
phases,  expected  at  the  beginning  or  just  prior  to  the  Lg  window,  are  not 
separated  or  distinct  on  the  records.  The  Pg  phases  are  most  strongly  excited 
for  the  Swiss  munitions  blast  and  Polish  rockburst,  and  Pn  is  strongest  for  the 
rockburst.  The  events  all  produce  quite  similar  maximum  amplitude  levels  at 
GRFO.  Fourier  spectral  analyses  were  performed  on  the  time  segments 
corresponding  to  each  of  the  identifiable  regional  phases  for  each  event,  as 
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Figure  13.  Locations  of  three  different  source  types  used  in  discrimination  analyses  for  station  GRFO 
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Figure  14.  Vertical-component  records  at  GRFO  from  three  different  source  types. 


well  as  for  noise  segments  preceding  the  Pn  signals  for  each.  Window 
segments  used  in  these  analyses  were  generally  12.5  seconds  in  duration 
except  in  one  or  two  cases  where  shorter  segments  were  used  to  avoid  some 
minor  data  glitches  which  contaminated  the  records  for  parts  of  the  regional 
phase  windows.  Figure  15  shows  signal  and  noise  spectra  for  each  of  the 
regional  phases  as  determined  from  the  GRFO  recording  of  the  8/15/81  Polish 
rockburst  and  is  typical  of  the  spectra  for  all  the  events.  The  S/N  levels  are 
generally  good  to  somewhere  in  the  6  -  8  Hz  range;  S/N  levels  are  lowest  for 
the  Pn  phases.  Even  though  the  S/N  is  much  greater  than  1.0,  earthquake 
spectra  tend  to  flatten  out  above  6  Hz  which  might  be  indicative  of  some  type  of 
noise  problem  at  higher  frequency. 

Previous  studies  of  regional  seismic  discrimination  (cf.  Bennett  et  al., 
1992)  have  demonstrated  the  potential  merits  of  comparing  relative  excitation  of 
different  regional  phases  for  use  in  identifying  source  types.  We,  therefore, 
decided  to  compute  the  spectral  ratios  of  the  individual  regional  phases  to  one 
another  -  i.e.  Lg/Pn  and  Lg/Pg.  We  show  here  the  Lg/Pg  ratios  because  of  the 
better  S/N  conditions  for  Pg  and  consequent  lower  danger  of  noise 
contamination  in  the  ratios.  These  ratios  for  the  three  events  recorded  at  GRFO 
are  plotted  in  Figure  16  for  the  band  0.5  -  6  Hz  in  which  S/N  is  generally  well 
above  one  for  the  individual  phases.  The  Lg/Pg  ratio  for  the  rockburst  tends  to 
be  highly  oscillatory  with  frequency  but  maintains  an  overall  level  just  below  a 
value  of  2  with  no  particular  evidence  of  increase  or  decrease  with  increasing 
frequency.  The  Lg/Pn  spectral  ratio  for  this  event  produced  quite  similar 
behavior  to  the  Lg/Pg  ratio.  Lg/Pg  for  the  earthquake  tends  to  be  less  oscillatory 
about  a  value  near  3  and  may  increase  slightly  toward  higher  frequencies.  The 
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Figure  15.  Example  of  signal  and  noise  spectra  for  the  regional  phases  recorded  on  the  vertical 
component  at  GRFO  from  the  8/15/81  Polish  rockburst. 
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Figure  16.  Lg/Pg  spectral  ratios  determined  at  GRFO  for  the  three  different  source  types. 


Lg/Pn  ratio  for  the  earthquake  was  somewhat  more  oscillatory  and  showed  a 
stronger  tendency  to  increase  with  frequency.  Lg/Pg  for  the  Swiss  blast  is  again 
highly  oscillatory  but  maintains  a  much  lower  level  near  0.7.  There  is  also  a 
tendency  in  the  ratio  to  decrease  with  frequency,  which  again  appears  to  show 
up  more  strongly  in  the  Lg/Pn  ratio  for  the  blast. 

In  Figures  17  and  18  we  compare  the  individual  phase  spectra,  Pg  and 
Lg,  observed  at  GRFO  for  the  events.  For  the  8/15/81  rockburst  and  the  4/14/92 
earthquake  in  Figure  17,  the  Pg  spectra  appear  generally  similar  in  level  and 
shape;  although  the  earthquake  shows  a  slight,  but  definite,  tendency  to  have 
more  higher  frequency  energy.  The  Lg  spectrum  for  the  rockburst  has  a  similar 
level  to  the  earthquake  in  the  1  -  2  Hz  band  but  drops  off  much  more  rapidly  at 
high  frequencies  than  the  corresponding  earthquake  spectrum.  In  Figure  18  the 
Pg  spectrum  for  the  Swiss  blast  is  strongly  enriched  (by  a  factor  of  about  3)  in 
the  band  1  -  2.5  Hz  and  slightly  enhanced  in  the  band  4.5  -  6  Hz  compared  to 
the  earthquake.  The  Lg  spectrum  for  the  Swiss  blast  and  the  earthquake  nearly 
coincide  up  to  about  3  Hz,  but  above  3  Hz  the  blast  Lg  spectrum  drops  off 
abruptly  to  a  level  3  or  more  below  the  earthquake.  Propagation  differences 
related  to  the  different  paths  into  GRFO  for  these  three  events  might  contribute 
to  the  variations  in  spectral  behavior,  but  source  mechanism  differences  might 
also  be  a  cause.  Additional  investigations  of  regional  propagation  effects  for 
this  area  and  theoretical  modelling  of  the  source  effects  could  help  resolve 
which  factors  are  significant. 

The  behavior  at  GRFO  of  the  P/P  spectral  ratios  between  different  event 
types  is  illustrated  in  Figure  19  for  Pg  phases.  The  rockburst  to  earthquake  ratio 
shows  a  clear  and  steady  decrease  with  frequency.  This  observation  appears 
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Figure  17.  Comparison  of  Pg  spectra  and  Lg  spectra  determined  at  GRFO  for  the 
rockburst  and  earthquake  from  Figure  13. 
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Figure  18.  Comparison  of  Pg  spectra  and  Lg  spectra  determined  at  GRFO  for  the 
explosion  and  earthquake  from  Figure  13. 


08/15/81  ROCKBURST  /  04/14/92  EARTHQUAKE 
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Figure  19.  Comparison  of  P/P  spectral  ratios  determined  at  GRFO  from 
the  different  event  types  from  Figure  1 3. 
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to  be  consistent  with  our  previous  observations  reported  above  for  teleseismic  P 
and  for  far-regional  P  from  German  and  some  Polish  rockbursts  relative  to 
earthquakes.  The  behavior  of  the  Pg/Pg  spectral  ratio  for  the  11/02/92 
explosion/4/14/92  earthquake  is  more  complex.  After  a  slow  decrease  with 
frequency  to  near  4  Hz,  there  appears  to  be  a  slight  increase  toward  higher 
frequencies.  The  Pn/Pn  spectral  ratios,  not  shown,  exhibit  a  steady  decrease 
with  frequency  over  the  entire  band  for  the  rockburst/earthquake  ratio;  while  the 
blast/earthquake  ratio  show  a  gradual  increase  with  frequency  over  the 
available  frequency  band.  These  latter  results  suggest  the  possibility  of  some 
discrimination  capability  between  different  source  types  in  the  regional  phase 
spectra,  but  additional  studies  are  needed  to  resolve  propagation  path 
influences  and  demonstrate  the  capability  for  larger  samples  of  events. 


48 


4.  South  African  Rockbursts 


4.1  Teleseismic  Signals 

As  noted  above,  many  large  rockbursts  occur  in  the  gold-mining  region  of 
South  Africa.  While  these  large  events  (Ml  -  5)  frequently  produce  good 
signals  at  teleseismic  stations  throughout  much  of  the  world,  there  is  little 
natural  earthquake  activity  or  man-made  explosions  of  similar  size  from  this 
source  region  which  can  be  used  to  provide  comparisons  for  use  in 
discrimination  studies.  Figure  20  shows  a  comparison  of  vertical-component 
records  at  three  teleseismic  AEDS  stations  for  a  large  South  African  rockburst 
(10/28/86,  5.2  mb)  and  the  largest  presumed  earthquake  (9/05/83,  4.7  mb)  from 
our  database.  As  can  be  seen  by  comparison  to  other  events  in  Table  2  above, 
this  earthquake  was  located  considerably  south  and  west  of  the  primary 
rockburst  zone.  There  also  appear  to  have  been  some  aftershocks,  including 
one  within  about  seven  seconds  after  the  main  event,  reported  in  the  NEIC 
catalog.  These  observations  tend  to  corroborate  the  identification  of  the  event 
as  an  earthquake. 

Comparing  the  signals  the  P  waveform  for  the  rockburst  appears  to  be 
significantly  more  complex  than  that  of  the  earthquake.  The  rockburst  P 
appears  to  include  several  oscillations  which  could  be  attributed  to  multiplicity 
in  the  source  or  possibly  reverberation  in  the  source  crust.  On  the  other  hand, 
the  earthquake  P  seems  to  be  a  rather  simple  pulse  (particularly  noticeable  at 
the  BC  station).  This  behavior  is  similar  to  that  described  in  Section  3. 1  above 
for  the  signals  from  Central  European  earthquakes  and  rockbursts  at 
teleseismic  stations.  The  second  pulse,  arriving  about  7  seconds  after  the  initial 
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Figure  20.  Comparison  of  vertical-component,  teleseismic  P  signals  at  three  AEDS  stations  from  a 
large  South  African  rockburst  and  a  nearby  presumed  earthquake. 


P  in  the  earthquake  record,  has  been  identified  as  P  from  an  aftershock, 
according  to  the  NEIC  interpretation.  If  the  event  is  indeed  an  aftershock,  it 
appears  that  the  mechanism  must  be  substantially  different  -  in  particular,  the 
pulse  polarity  seems  to  be  reversed  at  station  BC  and  possibly  at  the  other 
stations.  The  limited  waveform  data  that  we  have  looked  at  here  might  support 
an  alternate  interpretation  of  the  second  pulse  as  a  depth  phase.  If  the  phase  is 
pP,  the  focal  depth  for  the  earthquake  corresponding  to  the  7-second  delay  is 
about  20  km.  We  do  not  know,  at  this  time,  whether  other  observations, 
including  regional  data,  would  permit  this  alternative  interpretation. 
Nevertheless,  following  either  interpretation  appears  to  support  the  conclusion 
that  the  main  event  was  an  earthquake. 

4.2  Spectral  Comparisons  of  Teleseismic  P 

In  addition  to  simple  waveform  comparisons,  we  have  performed  spectral 
analyses  of  the  teleseismic  P  signals  recorded  at  the  AEDS  stations  for  several 
South  African  events.  These  analyses  were  similar  to  those  performed  on  the 
AEDS  teleseismic  signals  for  Central  European  events  as  described  in  Section 

3.2  above.  Many  of  the  teleseismic  AEDS  stations  are  at  very  large  epicentral 
distances  from  South  African  events  (65°  <  A  <  145°).  As  a  result,  many  of  the 
P  signals  from  South  African  events  have  experienced  greater  attenuation,  and 
the  useful  frequency  band  of  the  spectra  is  more  limited.  In  addition,  while  the 
signal  at  the  nearer  is  the  normal  teleseismic  P,  the  signal  at  the  more  distant 
stations  is  PKP.  Since  our  analyses  utilize  spectral  ratios  between  nearby 
events  at  the  same  station,  whether  we  are  using  P/P  or  PKP/PKP  ratios  should 
not  be  significant  to  our  interpretation.  We  have  analyzed  the  signals  from  three 
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of  the  largest  South  African  rockbursts  in  our  database:  1/25/89,  5.5  mb; 
9/26/90,  5.4  mb;  and  1 1/03/91,  5.1  mb-  Nine  common  stations  were  selected  for 
the  first  two  events,  and  six  for  the  latter  two  events. 

Fourier  spectra  for  7.5-second  P-wave  segments  and  similar  noise 
segments  were  calculated,  and  signal  versus  noise  analyzed  as  a  function  of 
frequency.  At  some  stations  signals  were  only  above  noise  out  to  about  3  Hz. 
P/P  spectral  ratios  were  computed  using  the  common-station  signal  spectra 
pairs  for  1/25/89/9/26/90  and  1 1  /03/9 1  / 9/26/90.  The  spectral  band  in  which 
S/N  ratios  were  generally  greater  than  one  was  again  used  to  limit  the  stations 
contributing  to  an  average  ratio  which  was  computed  for  each  event  pair.  The 
average  ratios  with  plus-  and  minus-one-sigma  bounds  about  the  mean  are 
shown  in  Figure  21.  In  both  cases  the  average  P/P  ratios  show  some 
fluctuations  with  frequency,  but  the  fluctuation  about  the  overall  trends  are  not 
too  great.  For  the  1/25/89/9/26/90  ratio,  the  ratio  oscillates  about  a  level 
somewhat  above  1.0.  There  also  appears  to  be  a  slight  decrease  in  the 
average  ratio  with  increasing  frequency  from  a  value  near  2  at  1  Hz  to  a  value 
near  1  at  4.5  Hz.  For  the  1 1  /03/9 1  / 9/2 6/9 0  P/P  spectral  ratio,  the  fluctuations 
are  somewhat  less  regular.  The  level  of  the  ratio  is  generally  below  one  over 
much  of  the  band  and  shows  a  tendency  to  increase  with  frequency  from  a 
value  near  0.5  at  1  Hz  to  a  value  near  1.0  at  3.5  Hz.  This  observation  of 
differences  in  the  frequency  content  of  teleseismic  P  for  rockbursts  at  nearby 
source  locations  appears  to  corroborate  our  preliminary  observations  (cf. 
Bennett  et  al.,  1993)  based  on  GDSN  station  measurements. 

Several  factors  may  contribute  to  the  observed  differences  in  P-wave 
spectra  and  P/P  spectral  ratios  even  though  the  source  types  are  similar.  The 
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stations  for  large  South  African  rockbursts. 


first  two  South  African  rockburst  events  are  in  roughly  the  same  location  (mine) 
while  the  third  is  almost  100  km  to  the  north.  It  is  unclear  to  what  extent  source 
location  differences,  such  as  different  site  response,  might  cause  the  observed 
variations  in  the  two  P/P  spectral  ratios.  In  addition,  the  differences  in  event 
sizes  is  not  great  but  might  be  responsible  also  for  some  frequency  and  level 
dependence  in  the  ratios.  The  smaller  event  would  be  expected  to  have  a 
somewhat  higher  corner  frequency,  so  that  the  P/P  spectral  ratio  of  a  smaller 
event  to  a  large  event  should  be  enhanced  at  higher  frequencies.  This,  in  fact, 
appears  to  be  consistent  with  our  observations  in  that  the  increasing  ratio  in 
Figure  21  corresponds  to  the  ratio  of  a  smaller  to  a  larger  event,  while  the 
decreasing  ratio  in  Figure  21  corresponds  to  the  ratio  of  a  larger  to  smaller 
event.  However,  other  factors  including  mechanism  variations  might  also 
contribute  to  the  observation.  For  example,  events  with  similar  body-wave 
magnitudes  could  be  expected  to  produce  different  P/P  spectral  ratios  because 
the  source  time  function  might  be  different  for  a  complex  roof  collapse  than  for  a 
simple  fracture  mine  tremor.  As  noted  above,  well-designed  theoretical  studies 
along  with  some  additional  analyses  of  observations  could  do  much  to  resolve 
some  of  the  outstanding  issues. 

With  regard  to  comparisons  of  teleseismic  P  spectra  between 
earthquakes  and  rockbursts  in  South  Africa,  we  have  used  the  P  signals 
recorded  at  the  three  AEDS  stations  from  the  two  events  described  above  in 
Section  4.1  (viz.  the  10/28/86  rockburst  and  the  9/05/83  earthquake.  The 
Fourier  spectra  were  computed  from  a  6-second  window  of  the  P-wave  (cf. 
Figure  20  above)  which  was  selected  to  exclude  the  P  signal  of  the  aftershock 
(or  pP  phase)  from  the  main  earthquake  P  segment.  We  again  computed  P/P 
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spectral  ratios  for  the  rockburst/earthquake  and  computed  an  average  ratio  over 
those  frequencies  in  which  the  S/N  level  was  determined  to  be  adequate  at 
each  station.  This  average  P/P  spectral  ratio  is  shown  in  Figure  22  over  the 
frequency  band  from  0.5  to  5  Hz.  The  ratio  is  fairly  consistent  between  stations 
with  each  showing  a  relative  maximum  in  the  vicinity  of  1  Hz  and  a  relative 
minimum  just  above  1.5  Hz.  The  ratio  fluctuates  somewhat  about  a  level  near 
3.0  at  lower  frequencies  consistent  with  the  0.5  magnitude  unit  difference 
between  the  two  events.  At  the  higher  frequency  end,  the  spectral  ratio  is  lower; 
but  we  are  less  confident  in  the  values  there  because  of  noise  contamination  at 
some  stations.  If  we  accept  the  lower  P/P  ratio  at  high  frequencies,  the 
observation  here  would  be  consistent  with  the  main  results  which  we  reported 
above  in  Section  III  for  the  Central  European  rockburst  and  earthquake 
comparisons.  However,  the  trends  here  are  not  strong,  and  we  would  like  to 
have  additional  corroboration. 

4.3  Regional  Observations 

Our  regional  signal  analyses  for  South  African  events  is  limited  to 
observations  from  the  near-regional  DWWSSN  station  SLR;  there  are  no  high- 
quality  stations  at  far-regional  distances  for  the  South  African  rockburst  zone. 
The  distance  to  SLR  from  events  in  the  South  African  rockburst  zone  is  from 
about  50  km  to  about  300  km,  and  most  regional  earthquakes  are  located  at 
somewhat  greater  distances.  For  a  database  including  31  mine  tremors  and  10 
regional  earthquakes  recorded  at  SLR,  Bennett  et  al.  (1993)  found  that  the  time- 
domain  measurements  of  maximum  P  versus  maximum  Lg  were  intermingled. 
This  is  illustrated  in  Figure  23,  extracted  from  that  previous  report.  The  Lg/P 
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Figure  22.  Average  P/P  spectral  ratio  determined  from  three  teleseismic  AEDS 
stations  for  a  large  South  African  rockburst  and  nearby  presumed  earthquake. 
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Figure  23.  Comparison  of  the  maximum  P  versus  maximum  Lg 
amplitudes  measured  from  time-domain  records  at  SLR  from  31  mine 
tremors  and  10  southern  Africa  earthquakes. 
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ratios  for  both  the  rockburst  and  earthquakes  from  South  Africa  have  values 
well-above  one  and  show  no  tendency  to  separate  relative  to  source  type.  We 
also  showed  in  the  previous  report  that  adjustments  for  propagation  path  length 
did  not  seem  to  make  much  difference  in  this  observation. 

In  the  discussion  of  the  teleseismic  P  observations  above  for  South 
Africa,  we  noted  apparent  variability  between  rockburst  events  from  essentially 
the  same  source  area.  Regional  observations  at  station  SLR  for  some  of  the 
smaller  rockburst  events  in  our  database  appear  to  corroborate  variability  in 
spectral  behavior  and  mechanism  between  nearby  mine  tremors.  Figure  24 
shows  a  comparison  of  the  Lg  spectra  for  two  mine  tremors  from  the  Orange 
Free  State  mining  area.  These  two  events  are  located  at  a  range  of  about  290 
km  from  SLR  and  are  separated  from  each  other  by  less  than  5  km. 
Considering  the  uncertainty  in  epicenter  location  within  the  region,  the  two 
sources  are  essentially  co-located.  The  Lg  spectra  for  the  two  events  have 
somewhat  similar  shapes.  The  difference  in  amplitude  level  in  the  frequency 
band  from  about  0.5  to  2  Hz  amounts  to  about  a  factor  of  three  which  is  just 
slightly  greater  than  the  0.4  magnitude  unit  difference  between  the  two  events 
(viz.  3.6  Ml  versus  3.2  Ml  for  the  1 1/19/81  and  5/24/83  rockbursts  respectively). 
The  Lg  spectrum  for  the  1 1/19/81  rockburst  decays  more  rapidly  than  that  for  the 
5/24/83  rockburst,  and  as  a  result  the  latter  event  has  relatively  more  high 
frequency.  This  may  be  partly,  but  probably  not  completely,  accounted  for  by 
the  shift  in  corner  frequency  toward  higher  frequencies  for  the  smaller  event. 

Another  observation  which  appears  to  support  variability  in  the 
mechanism  of  South  African  rockbursts  is  the  regional  P-wave  first  motions.  In 
Figure  25  we  show  the  expanded  first  motion  plots  at  station  SLR  for  five 
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SPECTRA  AT  SLR  FOR  EVENTS  5/24/83  AND  11/19/81 


59 


Orange  Free  State  mining  area  of  South  Africa. 
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Figure  25.  Initial  P  waveforms  at  station  SLR  for  mine  tremors  from  the  Far  West 
Rand  mining  district  of  South  Africa. 
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rockburst  events  from  the  Far  West  Rand  mining  district  (R  «  100  km).  In  two 
cases  the  first  motion  is  clearly  dilatational;  in  two  cases  it  is  clearly 
compressional;  and  one  case  is  ambiguous.  Mine  tremor  events  separated  by 
less  than  about  5  km  show  opposite  first  motions.  It  seems  most  likely  that  this 
observation  can  be  best  explained  by  a  change  in  source  mechanism  between 
the  events.  It  is  not  unreasonable  to  expect  that  such  a  source  mechanism 
difference  could  affect  the  spectral  behavior  of  the  transmitted  seismic  signals  in 
addition  to  the  first  motions. 

A  reasonable  interpretation  of  the  regional  observations  at  SLR  from 
events  in  the  South  African  rockburst  zone  is  that  there  is  some  degree  of 
variability  in  source  mechanism  between  events  from  a  rather  limited  source 
area.  The  mechanism  differences  are  the  most  likely  cause  of  the  observed  P- 
wave  first  motion  differences  and  may  explain  some  of  the  spectral  differences 
observed  in  regional  and  teleseismic  signals. 
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5.  Rockburst  Mechanisms  and  Other 
Discriminant  Measures 

5.1  Investigation  of  the  Ms  versus  mb  Discriminant 

A  teleseismic  discriminant  which  has  proven  to  be  extremely  valuable  for 
distinguishing  underground  nuclear  explosions  from  earthquakes  involves 
comparison  of  the  relative  excitation  of  long-period  (T  «  20  seconds)  surface 
waves  versus  short-period  (T  =  1  second)  teleseismic  body  waves  for  different 
events.  The  discriminant  measure  is  usually  formulated  in  terms  of  a 
comparison  of  surface-wave  magnitude,  Ms,  versus  body-wave  magnitude,  mb, 
and  has  been  established  as  one  of  the  most  robust  identification  techniques 
available  for  shallow  events  (cf.  OTA,  1992).  Although  this  discriminant  has 
been  most  often  utilized  over  the  years  for  large  teleseisms,  investigators  (cf. 
Marshall  and  Basham,  1972;  Taylor  et  al.,  1986)  have  also  suggested  its 
extension  to  smaller  events  at  regional  distances.  However,  the  value  of  the 
Ms-versus-mb  discriminant  for  application  to  small  events,  including  chemical 
blasts  and  rockbursts,  remains  generally  unestablished  for  the  variety  of 
potential  source  environments  of  interest  under  a  CTBT.  In  particular, 
calibration  studies  in  various  source  areas  and  tectonic  regions  would  be  useful 
in  determining  discriminant  characteristics  and  validity. 

In  the  current  investigation  we  have  applied  the  Ms/mb  discriminant  to 
rather  limited  event  samples  from  Central  Europe  and  South  Africa.  The 
samples  include  rockbursts  and  earthquakes  but  are  limited  by  the  fact  that  Ms 
values  are  seldom  reported  for  smaller  events.  Since  the  majority  of 
earthquakes  and  rockbursts  in  these  two  source  regions  are  small,  Ms 
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estimates  are  not  often  available.  This  appears  to  be  a  special  problem  for 
rockbursts  as  events  with  mb  or  Ml  greater  than  4.5  often  do  not  have 
associated  Ms  estimates. 

Figure  26  shows  a  sample  of  12  events  gleaned  from  the  NEIC  catalogs 
for  Central  Europe  and  South  Africa  for  which  Ms  measurements  were 
reported.  The  sample  includes  three  earthquakes,  and  the  remainder  are 
1  rockbursts.  In  the  figure  we  have  plotted  the  NEIC  reported  Ms  versus  NEIC  mb 
for  each  event.  For  reference,  we  show  lines  for  which  Ms  =  mb  and  for  which 
Ms  =  mb  -  1.0.  The  measurements  show  more  than  an  order  of  magnitude  of 
scatter.  With  regard  to  discrimination,  the  traditional  thinking  is  that  the  more 
that  event  Ms  values  fall  below  mb  values  the  more  likely  the  event  is  to  be  an 
explosion.  For  four  out  of  the  12  events  shown  in  Figure  26,  the  Ms  values  are 
within  0.3  magnitude  units  of  mb;  they  include  two  of  the  earthquakes. 
However,  the  remaining  eight  events  have  Ms  values  more  than  one-half 
magnitude  unit  lower  than  mb  and  include  four  events  with  Ms  values  more 
than  one  magnitude  unit  lower  than  mb- 

Figure  27  shows  a  similar  plot  of  magnitude  measurements  derived  by 
averaging  AEDS  station  magnitudes  for  a  sample  of  ten  events.  It  should  be 
noted  that  some  of  the  same  events  as  in  the  preceding  figure  are  included  here 
and  that  no  station  corrections  were  applied  in  computing  the  event  averages. 
In  this  case  the  Central  European  earthquake  (4/13/92)  falls  near  the  Ms  =  mb 
line,  but  all  of  the  rockbursts  appear  to  be  more  explosion-like.  Comparing  just 
the  Central  European  rockburst  (3/13/89)  and  the  Central  European 
earthquake,  the  earthquake  appears  quite  normal  (Ms  *  mb);  but  the  rockburst 
looks  very  explosion-like  (Ms  <  mb  -1.0).  Considering  that  both  of  the  latter 
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Figure  26.  M$  versus  mb  measurements  reported  in  NEIC  catalog  for 
selected  rockbursts  and  earthquakes  from  Central  Europe  and  South 
Africa. 


Figure  27.  Ms  versus  mb  determined  from  AEDS  station 
measurements  for  selected  rockbursts  and  an  earthquake  from 
Central  Europe  and  South  Africa. 


events  come  from  essentially  the  same  source  region,  it  seems  unlikely  that 
path  corrections  would  do  much  to  compensate  this  difference.  The  observation 
suggests  that  rockbursts  may  tend  to  look  more  explosion-like  than  earthquakes 
with  respect  to  the  Ms/mb  discriminant.  However,  this  conclusion  is  rather 
weak  because  of  the  limited  data  sample  from  which  it  is  drawn,  and  it  should 
probably  be  tested  for  a  larger  sample. 

For  the  South  African  events,  we  did  try  to  assemble  a  somewhat  larger 
sample  of  Ms  measurements  using  data  recorded  at  station  BCAO  (A  «  32°).  In 
our  previous  report  (Bennett  et  al.,  1993)  we  showed  evidence  of  fairly  strong 
long-period  Rayleigh  waves  at  BCAO  from  several  events  in  our  South  Africa 
rockburst  database.  For  a  sample  including  18  rockburst  events  with  4.2  <  mb  < 
5.2,  we  computed  surface-wave  magnitudes  from  the  long-period  Rayleigh 
waves  recorded  at  BCAO.  The  waveform  segments  were  first  filtered  in  a 
passband  from  0.04  Hz  to  0.06  Hz  to  isolate  signals  with  periods  in  the  vicinity 
of  20  seconds.  Peak-to-trough  vertical-component  ground  motion  amplitudes 
determined  from  the  maximum  amplitude  in  the  Rayleigh  wave  group  velocity 
window  were  then  used  in  the  Prague  formula  (Bath,  1979)  to  calculate  Ms: 

Ms  =  log  (A/T)  +  1 .66  log  A  (4) 

where  A  is  the  measured  amplitude  in  nanometers,  T  is  period  in  seconds,  and 
A  is  epicentral  distance  in  degrees.  The  results  were  plotted  in  Figure  28 
versus  the  NE1C  reported  body-wave  magnitude,  mb-  For  two  events  the 
Rayleigh  waves  appeared  to  have  some  interference  from  other  events;  these 
are  marked  with  downward  arrows  suggesting  the  corresponding  Ms  estimate 
may  be  an  upper  bound.  Although  the  Ms-versus-mb  measurements  in  Figure 
28  show  considerable  scatter,  they  appear  to  be  consistently  explosion-like  with 
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Figure  28.  Mg  determined  from  long-period  Rayleigh  wave 
measurements  at  station  BCAO  versus  m  b  reported  by  NEIC  for 
selected  South  African  rockbursts. 
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Ms  values  in  all  cases  more  than  one  magnitude  unit  below  mb  and,  in  some 
cases,  as  much  as  two  magnitude  units  below  mb.  It  should  be  noted  here, 
however,  that  measurements  are  at  a  single  station,  and  no  station  corrections 
have  been  applied.  If  the  transmission  path  into  BCAO  is  poor  or  the  station  lies 
on  a  node  of  the  Rayleigh-wave  radiation  pattern,  low  Ms  values  might  be 
expected.  Additional  investigations  using  other  stations  for  some  of  the  larger 
events  and  predictions  based  on  knowledge  of  the  focal  mechanisms  could 
help  to  resolve  the  cause  for  the  low  Ms  measured  for  these  South  African 
rockbursts. 

5.2  Rockburst  Mechanisms 

In  our  previous  report  (Bennett  et  al.,  1993)  we  described  some  of  the 
mechanisms  associated  with  rockburst  activity.  Six  distinct  models  for  mining- 
induced  seismic  events  were  defined  by  Hasegawa  et  al.  (1989).  These 
models  are  (1)  cavity  collapse,  (2)  pillar  burst,  (3)  tensile  failure  in  the  cap  rock, 
(4)  normal  faulting  above  the  stope  face,  (5)  thrusting  above  or  below  the 
excavation,  and  (6)  shallow,  horizontal  thrusting  along  layering  planes  above 
the  mine  roof.  These  various  models  would  be  expected  to  differ  with  respect  to 
energy  release  and  the  mechanisms  of  seismic  wave  generation.  Individual 
rockfalls  or  pillar-bursts  would  be  expected  to  normally  release  only  small 
amounts  of  seismic  energy  and  to  have  simple  vertical-force  or  dipole 
mechanisms  unlike  the  double-couple  mechanisms  experienced  in  some  other 
rockbursts  and  typical  earthquakes.  Although  several  of  these  models  may  act 
in  combination  to  produce  quite  complex  seismic  mechanisms,  for  many 
rockbursts  in  many  different  source  regions  the  mechanisms  appear  to  be 
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dominated  by  shear  failure  similar  to  that  in  earthquakes.  Furthermore,  in  many 
mining  areas  the  ambient  tectonic  stress  appears  to  contribute  to  the 
mechanism,  enhancing  the  source  energy  release  and  controlling  the  rockburst 
process.  As  a  result,  many  rockbursts  throughout  the  world  would  be  expected 
to  behave  in  most  aspects  similar  to  shallow  earthquakes  with  regards  to 
seismic  discrimination  characteristics.  However,  unusual  events  with  more 
complex  mechanisms  sometimes  do  occur  and  may  be  found  in  essentially  the 
same  source  region  where  the  simpler  mechanism  dominates.  Some  of  these 
complex  events  have  been  large,  as  will  be  described  below;  and  there  is 
sometimes  a  close  tie  to  mining  practice. 

The  mechanisms  of  rockbursts  in  South  Africa  have  been  studied  and 
analyzed  more  than  those  in  any  other  area.  The  attention  they  have  received 
can  be  attributed  to  several  factors  including  the  size  and  frequency  of  the  mine 
tremors  and  the  severe  consequences  of  some  events.  In  addition,  the 
substantial  seismic  monitoring  effort  in  the  vicinity  of  South  African  mines 
provides  excellent  capability  to  resolve  mechanisms  and  analyze  the  influences 
of  mining  practice.  One  result  of  this  is  that  there  is  some  tendency  for 
experience  gained  from  South  African  mine  tremors  to  dominate  understanding 
of  the  rockburst  problem  in  other  parts  of  the  world.  In  recent  years  more 
detailed  analyses  have  been  made  of  rockburst  mechanisms  from  mining  areas 
in  other  geographical  regions  (e.g.  Central  Europe  and  Canada)  which  also 
have  extensive  records  of  such  activity.  In  many  cases  these  latter  studies 
confirm  the  South  African  experience,  but  distinctive  features  also  have  been 
found  in  some  cases. 
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In  South  Africa  observations  from  large  mine-induced  events  support  the 

conclusion  that  such  events  are  primarily  associated  with  shear  failure  on  faults 

or  zones  of  weakness  in  the  vicinity  of  the  excavation  (cf.  McGarr,  1971; 
Spottiswoode  and  McGarr,  1975;  McGarr  et  al.,  1979).  These  large  events  are 
usually  closely  associated  with  the  prevalent  tectonic  stress  conditions  in  the 
region  surrounding  the  mine,  so  that  the  seismic  mechanisms  represent  the 
response  of  fractures  or  weaknesses  adjacent  to  the  excavation  which  are 
favorably  oriented  in  the  ambient  stress  field.  Gay  et  al.  (1984)  argue  that 
mining-induced  energy  cannot  possibly  account  for  a  significant  part  of  the  total 
energy  released  during  a  large  event.  They  conclude  that  the  large  mine 
tremors  in  South  Africa  draw  from  residual  strain  energy  related  to  the  regional 
tectonics. 

Many  small  events,  with  magnitudes  in  the  range  1  <  Ml  £  3  also 
produce  a  quadrantal  distribution  of  first  motions,  as  measured  on  near-source 
arrays,  which  suggest  that  the  events  are  the  result  of  shear  failure  (cf. 
Spottiswoode,  1984;  Potgieter  and  Roering,  1984).  However,  there  have  been 
some  exceptions  to  this  behavior  even  in  South  Africa.  Rorke  and  Roering 
(1984)  found  that  mine-tremor  events  with  magnitudes  in  the  range  -1.0  <  Ml  £ 
2.8  at  the  Doornfontein  Gold  Mine  did  not  produce  quadrantal  first-motion 
patterns  but  had  only  two  lobes.  In  most  events  occurring  close  to  the  mine 
stope  face,  a  rarefactional  lobe  occurred  on  the  stope  side  of  the  event  and  a 
compressional  lobe  on  the  unmined  side.  Exceptions  to  the  double-couple 
mechanism  for  rockbursts  have  also  been  found  in  Central  Europe  (cf.  Sileny, 
1989;  Gibowicz  et  al.,  1990)  and  in  the  western  U.S.  (cf.  Wong  et  al.,  1989; 
Taylor,  1992).  In  fact,  the  large  German  rockburst  which  will  be  described  more 
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completely  later  in  this  section  also  may  be  an  exception  to  the  simple  shear 
mechanism. 

In  our  previous  report  (Bennett  et  alM  1993)  we  described  the  more 
typical  rockburst  model  as  a  seismic  event  triggered  when  the  stresses  are 
redistributed  on  a  tightly  confined  fault  usually  in  the  vicinity  of  the  excavation. 
This  stress  redistribution,  which  is  caused  by  the  mining  activity,  induces  a 
frictional  instability  on  the  fault  causing  a  shear  failure  to  spread  on  the  fault 
plane  and  radiating  seismic  waves.  Such  a  mechanism  is  equivalent  to  a 
double  couple  and  produces  the  typical  quadrantal  radiation  pattern.  The 
displacement  spectrum  from  this  type  of  source  consists  of  a  flat  low-frequency 
trend,  a  descending  intermediate  trend  with  slope  proportional  to  co’2  to  co'3  and 
an  upper  frequency  limit  controlled  by  either  attenuation  or  source  properties 
(cf.  Madariaga,  1979).  Studies  of  the  spectral  behavior  of  South  African  mine 
tremors  by  McGarr  et  al.  (1990)  and  Battis  (1992)  suggest  that  this  double¬ 
couple  rupture  model  satisfactorily  explains  the  events  which  they  analyzed. 
However,  in  some  cases  there  were  indications  in  their  results  that  the  rupture 
process  was  more  complex  probably  because  of  multiple  rupture  events. 

In  a  more  general  approach  to  the  problem  of  mining-induced  seismic 
events,  Kuhnt  et  al.  (1989)  distinguished  two  different  types  of  rockbursts:  (1)  a 
tectonic  type  related  to  induced  fault  movement  in  the  ambient  tectonic  stress 
field  (similar  to  the  typical  South  African  mechanism  just  described),  and  (2)  a 
static  type  which  is  more  closely  connected  to  mining  activity  and  not  dependent 
on  tectonic  stress.  For  this  latter  model  Kuhnt  et  al.  formulated  representations 
of  the  process  in  terms  of  both  barrier  and  asperity  models.  In  the  barrier  model 
the  surroundings  of  the  pillars  act  as  barriers  to  propagation  of  a  crack  outward 
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from  the  stope  area.  In  the  asperity  mode!  the  pillars  form  asperities  which 
concentrate  stresses  and  which  ultimately  fail  under  the  induced  stress  to 
produce  the  tremor.  In  either  model,  multiple  pillars  may  be  affected  in  a  single 
rockburst  and,  depending  on  their  locations  and  properties,  the  resulting 
rockburst  process  may  be  quite  complex.  The  studies  indicate  that  the  mining 
practice  can  to  some  degree  control  the  process  and  lead  to  variations  in  the 
source  between  events.  It  has  not  been  determined  to  what  extent  the 
complexity  in  the  model  may  contribute  to  the  observed  seismic  signals.  Some 
of  the  studies  already  cited  suggest  that  these  kinds  of  complexity  are  most 
likely  to  affect  small  events.  However,  there  may  be  some  mining  situations 
where  large  events  may  have  added  complexity. 

5.3  Mechanism  of  the  1989  Volkershausen,  Germany 
Rockburst 

In  several  preceding  sections  of  this  report,  we  showed  observations 
from  a  large  rockburst  in  Central  Europe.  This  event  occurred  on  March  13, 
1989  near  Volkershausen,  Germany  and  had  a  magnitude  of  5.4  mb  and  focal 
depth  near  1  km.  The  rockburst  occurred  at  the  Ernst  Thaelmann  Potash  Mine 
and  produced  considerable  damage  to  buildings  in  the  area  and  some  loss  of 
life  to  miners,  this  event  was  widely  recorded  on  local  and  regional  seismic 
networks  as  well  as  at  far-regional  and  teleseismic  stations.  A  few  reports  have 
been  published  regarding  the  mining  conditions  and  seismic  measurements 
associated  with  this  event  (e.g.  Ahorner,  1989;  Bormann,  1992;  Bowers  and 
Pearce,  1992;  Minkley,  1993).  We  review  here  the  source  models  and  some 
apparently  conflicting  information  observed  from  the  observed  seismic  signals. 
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Unlike  the  typical  South  African  mine  tremor,  the  Volkershausen 
rockburst  is  thought  to  have  been  produced  by  the  nearly  simultaneous 
collapse  of  multiple  mine  pillars  (cf.  Ahorner,  1989;  Minkley,  1993).  This  event 
was  apparently  triggered  when  an  attempt  was  made  to  mine-out  some  remnant 
pillars  in  the  mine  without  giving  appropriate  engineering  attention  to  effects  on 
the  surrounding  stress  field.  However,  the  consequence  was  far  greater  than 
failure  of  a  few  of  the  surrounding  pillars.  In  this  case  it  has  been  reported  that 
about  3200  pillars  collapsed  within  a  very  short  time  span  (cf.  Bormann,  1992). 
Figure  29  shows  a  schematic  model  of  the  Volkershausen  rockburst  along  with 
an  approximate  energy  calculation  to  estimate  the  energy  released  in  the  event. 
The  pillar  collapse  occurred  at  depths  of  about  850  -  900  m.  Minkley  (1993) 
shows  results  obtained  from  displacement  meters  at  the  surface  above  the  mine 
which  indicate  that  a  large  area  subsided.  The  subsidence  area  appears  to 
extend  over  the  surface  for  2  -  3  km  laterally.  The  displacement  meter  results 
further  indicate  that  maximum  subsidences  at  the  earth's  surface  were  from  0.8 
- 1.0  m  and  that  the  process  was  generally  completed  within  a  few  seconds. 

In  Figure  29  we  performed  a  simplified  calculation  to  determine  the 
potential  energy  loss  due  to  subsidence  by  0.7  m  of  a  block  approximately  3  km 
by  3  km.  We  also  calculated  the  seismic  energy  from  the  Gutenberg-Richter 
relationship  associated  with  a  seismic  event  having  a  body-wave  magnitude  of 
5.4  mb-  The  calculation  suggests  that  the  potential  energy  available  from  the 
mine  collapse  process,  in  this  case,  is  quite  large.  A  seismic  efficiency  of  only 
about  0.5  percent  would  be  needed  to  produce  the  observed  body-wave 
magnitude  for  this  event.  Even  if  the  average  subsidence  was  assumed  to  be 
only  half  as  great  or  the  lateral  extent  was  reduced  by  50  percent,  reasonable 
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Figure  29.  Schematic  of  the  1989  Volkershausen,  Germany  rockburst  (  5.4  m  t>)  mechanism  and  the 
associated  energy  relations. 


seismic  efficiencies  in  the  1  -  2  percent  range  could  still  be  obtained.  This 
calculation  indicates  that,  in  this  case,  the  mine  collapse  was  probably  large 
enough  to  produce  the  observed  magnitude  of  5.4  mb  without  having  to  draw 
on  tectonic  strain  energy.  Similar  calculations  of  this  type  could  be  performed 
easily  to  estimate  the  extent  of  mine  collapse  required  to  produce  events  of 
various  magnitudes.  More  sophisticated  model  studies  would  also  be  useful  for 
analyzing  details  of  this  and  other  rockburst  processes  and  to  estimate  the 
expected  ground  motions  at  local  and  regional  distances  from  such  events. 

With  regard  to  seismic  signal  generation,  Ahorner  (1989)  has  suggested 
that  the  Volkershausen  rockburst  source  could  be  represented  as  an  implosion. 
For  the  collapse  of  individual  pillars,  the  vertical  motion  vectors  above  and 
below  are  directed  inward  toward  the  pillar.  Horizontal  motions  are  directed 
outward  from  the  pillar,  but  these  motions  are  into  the  open  rooms  between 
pillars  and  do  not  transmit  as  seismic  waves.  The  model  is  not  strictly  an 
implosion  (for  which  horizontal  motions  would  also  be  directed  inward),  but 
rather  an  inward  vertical  dipole.  The  P— wave  first  motions  produced  by  this 
source  model  would  be  expected  to  correspond  to  dilatations  at  all  stations.  To 
test  the  model  Ahorner  studied  the  P-wave  signals  at  a  regional  network  of 
about  50  seismic  stations  out  to  ranges  of  380  km.  The  sense  of  the  first 
motions  observed  at  these  stations  are  shown  on  the  left  of  Figure  30,  as  plotted 
by  Ahorner.  All  the  first  motions,  corresponding  to  Pg  and  Pn  phases,  are 
dilatational  at  all  azimuths  which  apparently  agrees  with  the  proposed 
"implosional"  source  model. 

As  a  further  test  of  the  seismic  mechanism  of  this  event,  we  analyzed  the 
first  motions  for  all  stations  reported  in  the  ISC  Bulletin.  These  observations 
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Figure  30.  Comparison  of  P-wave  first  motions  for  the  March  13, 1 989  Volkershausen  rockburst  reported  by 
Ahorner  (1989)  from  observations  at  a  near-regional/regional  station  network  (left)  and  as  reported  in  the  ISC 
Bulletin  from  all  stations  (right).  Numerous  compressional  first  motions  at  ISC  stations  may  be  indicative  of 
signal-to-noise  limitations  which  prevent  detection  of  the  initial  pillar  failure. 


include  114  stations  at  ranges  from  A  <  1°  to  A  >  83°  and  incorporate  some,  but 
not  all,  of  the  stations  used  by  Ahorner.  These  first  motions  are  plotted  on  the 
right  of  Figure  30.  As  can  be  seen  in  the  plot,  both  compressional  and 
dilatational  first  motions  have  been  reported  from  these  stations.  Of  the 
reported  first  motions  62  are  dilatational  and  52  are  compressional.  Only  7  of 
the  35  Pn  phases  (reported  for  A  <  8°)  are  compressional,  and  only  7  of  the  30 
P  phases  reported  for  A  >  30°  are  dilatational.  There  appears  to  be  no  clear 
mechanism  pattern  in  the  distribution  of  the  first  motions.  One  explanation  for 
the  discrepancy  between  the  left  and  right  plots  in  Figure  30  could  be 
inaccurate  first  motion  reports,  presumably  from  the  ISC  Bulletin.  Such 
problems  arise  due  to  inconsistent  identifications  of  the  first  motions  in  low  S/N 
conditions.  However,  we  cannot  asses  the  validity  of  this  argument  without 
additional  knowledge  of  the  actual  waveforms  which  is  beyond  the  scope  of  the 
current  investigation.  A  related  explanation  of  the  distribution  of  first  motions  in 
the  right-hand  plot  might  be  that  the  Volkershausen  rockburst  was  a  complex 
event.  Near-in  stations  and  those  with  high  S/N  conditions  would  have  reported 
the  initial  “implosional"  phase  of  the  collapse.  However,  at  some  slightly  later 
time,  the  more  distant  teleseism ic  stations  and  those  with  lower  S/N  conditions 
might  have  seen  a  second,  possibly  larger,  event.  This  second  event  could 
have  been  associated  with  an  outward  vertical  dipole  phase  caused  by  closure 
between  the  roof  and  floor  of  the  mine  or  by  some  tectonic  stress  release.  The 
latter  might  also  explain  the  normal,  dip-slip  double-couple  source  reported  by 
Bowers  and  Pearce  (1992)  for  this  event  based  on  relative  amplitude  analyses 
of  initial  and  secondary  phases  at  two  teleseismic  station  arrays.  The  additional 
theoretical  studies  of  the  mechanism  of  this  rockburst,  as  noted  above,  along 
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with  detailed  analyses  of  the  observed  regional  and  teleseismic  signals  would 
be  useful  in  resolving  the  processes  contributing  to  the  Volkershausen 
rockburst. 

5.4  Implications  for  Identification  and  Control 

In  our  previous  report  (Bennett  et  al.,  1993)  we  noted  that  many 
rockbursts  have  mechanisms  similar  to  earthquakes  and  that,  in  general, 
discrimination  procedures  which  are  effective  for  earthquakes  should  also  work 
for  rockbursts.  It  was  further  suggested  that,  since  many  rockbursts  are  small, 
regional  discrimination  techniques  may  be  key  to  their  identification.  We 
considered  there  a  variety  of  traditional  discriminants  including  focal  depth,  Ms 
versus  mb,  sense  of  first  motion,  spectral  scalloping,  and  S  or  Lg  versus  P 
excitation.  Additional  comments  on  the  conclusions  presented  there  are 
appropriate  in  the  light  of  the  findings  of  this  current  report. 

First,  with  regard  to  focal  depth,  rockbursts  and  underground  nuclear 
explosions  occur  at  about  the  same  focal  depths  while  earthquakes  are  often 
located  much  deeper.  Depth  estimation  techniques,  such  as  depth  phases,  P- 
wave  location,  and  the  Wadati  technique  using  S  -  P  times,  would  not  be 
expected  to  be  able  to  reliably  determine  any  significant  differences  in  focal 
depths  between  rockbursts  and  underground  nuclear  explosion  tests.  So,  in 
this  respect  rockbursts  behave  differently  from  earthquakes  and  may  complicate 
utilization  of  some  standard  nuclear  explosion  discrimination  techniques  which 
rely  on  focal  depth  estimation. 

Second,  with  regard  to  the  Ms  versus  mb  discriminant,  we  noted  above 
that  Ms  is  not  often  determined  for  rockburst  events.  Many  of  the  events  are 
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small  and  not  very  well  recorded  on  long-period  instruments.  This  might  be  due 
to  a  reduced  excitation  of  long-period  (20-second)  Rayleigh  waves  for  the  small 
events.  For  even  the  largest  rockburst  events,  we  found  that  the  long-period 
Rayleigh  waves  were  relatively  weak.  Thus,  Ms  measurements  reported  for 
Central  European  and  South  African  rockbursts  were  frequently  on  the  order  of 
one-magnitude  unit  lower  than  mb-  Furthermore,  the  Ms  values  which  we 
estimated  at  BCAO  for  our  selected  sample  of  South  African  rockbursts  were 
from  one  to  two  magnitude  units  less  than  mb-  The  rockbursts  then  appear  to 
be  explosion-like  with  respect  to  the  Ms  versus  mb  discriminant  although  the 
degree  of  the  Ms  -  mb  difference  might  be  exaggerated  by  the  small  number  of 
stations  and  lack  of  station  corrections  in  our  magnitude  estimates.  In  any  case, 
application  of  the  Ms/mb  discriminant  to  small  events  including  rockbursts 
would  appear  to  require  some  assessment  of  these  excitation  differences  and 
enhancement  of  Ms  measurement  capability  probably  using  a  network  of 

regional  stations. 

Two  other  discrimination  methods,  first  motion  and  complexity,  appear  to 
offer  some  promise  for  distinguishing  rockbursts  from  other  source  types.  P- 
wave  first  motions  from  rockbursts  would  be  expected  to  differ  from  the  uniform 
compressions  expected  from  explosions.  Shear  rockbursts  should  produce 
quadrantal  radiation  patterns  similar  to  earthquakes;  however,  other  rockburst 
sources  may  be  less  regular.  In  particular,  we  saw  examples  of  rockbursts  with 
very  complex  first  motion  patterns  which  would  not  fit  an  earthquake  rupture 
model.  Furthermore,  P  waves  from  rockbursts  frequently  appear  to  be  emergent 
and  complex,  which  makes  first  motion  sometimes  hard  to  determine.  The 
complexity  of  the  rockburst  signal  does,  however,  lend  itself  to  an  alternate 
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discriminant  measure.  More  complex  signals  in  earthquakes  and  rockbursts 

should  provide  some  distinction  from  relatively  simple  nuclear  explosion 
sources.  In  fact,  we  found  evidence  in  this  study  of  complexity  in  rockburst 
signals  even  greater  than  in  earthquakes.  However,  additional  testing  of 
proposed  parametric  measures  of  complexity  are  needed  to  verify  this 
discriminant  for  rockbursts  as  well  as  earthquakes. 

With  regard  to  regional  signals,  we  found  that  S/P  spectral  ratios  fro 
rockbursts  tended  to  be  lower  than  those  of  earthquakes  but  not  as  low  as  those 
of  blasts.  In  addition,  there  was  evidence  of  differences  in  the  Lg/P  ratio 
behavior  as  a  function  of  frequency  for  the  different  source  types.  Earthquake 
and  rockburst  ratios  tended  to  increase  with  frequency  while  blast  ratios 
decreased.  We  also  found  that  regional  P  for  the  blasts  appeared  to  be 
enhanced  at  high  frequencies  compared  to  other  source  types.  Additional 
studies  are  needed  to  confirm  this  behavior  in  different  source  areas  and  for 
larger  event  samples;  however,  these  initial  results  offer  promise  that  regional 
spectra  may  provide  a  valuable  discriminant  measure  for  different  source  types. 

Finally,  a  particularly  interesting  aspect  of  rockbursts  bearing  on  the 
nuclear  explosion  test  evasion  problem  is  that  rockburst  sources  appear  to  be 
controllable.  The  idea  of  using  a  natural  earthquake  to  conceal  a  clandestine 
nuclear  explosion  test  was  originally  dismissed  because  the  potential  evader 
has  no  control  over  the  timing,  location,  and  size  of  the  earthquake  and, 
therefore,  cannot  plan  for  his  test.  However,  there  is  accumulating  evidence  in 
the  mining  engineering  literature  that  rockburst  sources  can  be  predicted  to 
some  degree  and  controlled.  In  our  previous  report  (Bennett  et  al.,  1993)  we 
noted  that  mining  engineers  have  developed  finite  element  techniques  which 
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enable  accurate  stress  analysis  and  prediction  in  the  rock  surrounding  mining 
excavations.  These  and  corroborative  monitoring  programs  provide  the 
capability  to  predict  rockburst  occurrence  and  level  with  some  accuracy  within 
specific  mining  areas.  Case  histories  have  been  reported  where  addition  or 
removal  of  supports  were  used  to  reduce  stress  level  below  some  critical 
threshold  or  to  increase  it  to  the  point  of  inducing  a  rockburst  in  some  selected 
magnitude  range. 

Therefore,  a  viable  evasion  scenario  appears  to  be  that  a  small  or 
decoupled  nuclear  explosion  test  might  be  concealed  in  a  triggered  rockburst 
event.  Given  the  widespread  nature  of  rockburst  occurrence  in  mines,  such  a 
scenario  might  be  carried  out  in  any  number  of  countries.  Assuming  that  even 
large  rockbursts  might  be  predictable  and  controllable,  the  companion  nuclear 
explosion  could  be  not  only  difficult  to  identify  but  its  seismic  signals  might  not 
even  be  detected.  Theoretical  model  studies  could  be  useful  in  evaluating  this 
scenario  and  identifying  possible  monitoring  strategies  for  specific  source 
areas. 


81 


6.  Summary  and  Conclusions 


6.1  Research  Summary 

As  a  result  of  changing  emphases  in  the  field  of  nuclear  explosion 
monitoring  in  recent  years,  alternate  types  of  seismic  sources  have  taken  on 
new  significance.  The  research  described  here  has  focussed  on  rockbursts,  or 
mining-induced  tremors,  because  they  represent  a  significant  class  of  events 
which  has  received  little  prior  attention  with  regards  to  seismic  discrimination. 
Rockbursts  are  important  to  the  seismic  discrimination  problem  because  they 
are  frequent,  occur  in  aseismic  as  well  as  seismic  areas,  may  be  large,  and  can 
be  controlled.  This  research  program  was  designed  to  investigate 
characteristics  of  rockbursts  which  might  help  facilitate  their  identification  in  the 
context  of  nuclear  explosion  test  monitoring. 

Although  rockbursts  occur  in  many  areas  throughout  the  world,  for  this 
study  we  focussed  on  two  particular  regions:  Central  Europe  and  South  Africa. 
In  both  of  these  source  regions,  there  have  been  frequent  reports  of  rockbursts 
and  several  of  these  have  been  large.  As  a  result,  good  waveform  data  from 
high  quality  stations  can  be  obtained  for  several  events  from  these  areas.  For 
both  Central  Europe  and  South  Africa,  we  attempted  to  assemble  a  database  of 
representative  events,  including  rockbursts  as  well  as  contrasting  source  types, 
which  could  be  used  in  discrimination  analyses.  Digital  seismic  waveforms 
were  collected  for  60  events  from  Central  Europe  and  79  events  from  South 
Africa.  The  waveforms  included  high-quality  data  from  SRO,  ASRO,  DWWSSN, 
ARPA  regional  arrays,  Soviet  IRIS,  GSETT-2,  and  AEDS  stations.  The  events 
covered  a  magnitude  range  from  2.0  to  5.5  mb-  In  all  cases  seismic  records 
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with  good  signal-to-noise  conditions  were  found  for  at  least  one  station.  The 
data  for  larger  events  included  many  teleseismic  waveforms  from  a  global 
network  of  stations,  while  the  smaller  events  provided  useful  data  at  only  a  few 
far-regional  or  regional  stations. 

For  these  data  a  variety  of  time-domain  and  spectral  analyses  were 
performed  on  the  waveforms  for  selected  events.  For  teleseismic  signals  we 
focussed  on  the  P-wave  segments  and  the  long-period  Rayleigh  waves  at 
selected  stations.  We  reviewed  the  complexity  of  the  P  waves  and  compared  P- 
wave  spectral  measurements  between  events.  We  also  analyzed  the  Ms* 
versus-mb  discriminant  for  larger  events  for  which  data  were  available.  For  the 
far-regional  and  regional  waveforms,  we  investigated  the  relative  excitation  of 
P,  S,  and  Lg  signals  in  both  the  time  and  spectral  domains  for  selected  events. 
We  also  looked  at  sense  of  first  motions  as  a  possible  discriminant  for  the 
rockbursts.  In  these  studies  we  attempted  to  select  events  and  stations  which 
would  reduce  to  some  degree  propagation  path  effects  and,  therefore,  provide 
indications  of  source  differences.  Comparisons  were  made  between  different 
source  types  to  identify  potential  distinguishing  characteristics  of  the  signals 
and  between  similar  source  types  to  identify  variability  in  the  source  or  possible 
path  influences. 

In  addition  to  the  data  analyses,  we  investigated  the  mechanisms  and 
control  of  rockburst  events.  In  these  investigations  we  reviewed  the  technical 
literature  to  identify  common  and  contrasting  features  in  the  source  mechanisms 
of  rockbursts.  We  then  sought  to  identify  what  implications  different 
mechanisms  and  source  characteristics  would  have  on  some  typical 
discriminant  measures.  As  part  of  this  review  of  rockburst  mechanisms,  we  also 


83 


considered  predictability  and  control  of  rockburst  events  by  mining  practice  and 
the  implications  of  such  capability  for  clandestine  nuclear  explosion  testing. 

6.2  Principal  Conclusions 

One  important  finding  of  this  research  program  is  that  rockbursts  occur  in 
hardrock  mining  areas  throughout  the  world.  In  many  areas  the  events  are 
frequent,  and  they  may  be  large.  The  widespread  nature  of  this  kind  of  activity 
indicates  that  they  could  pose  a  significant  problem  for  CTBT  or  proliferation 
monitoring.  In  fact,  in  many  parts  of  the  world  mining-induced  events  may 
present  more  of  a  problem  to  discrimination  monitoring  than  commercial 
blasting  operations. 

Among  the  techniques  for  distinguishing  rockbursts  from  other  source 
types,  the  spectral  differences  in  seismic  signals  appear  to  be  most  promising. 
The  Lg/Pg  spectral  ratios  observed  at  regional  station  GRFO  for  different  Central 
European  source  types  at  similar  regional  distances  suggest  that  the  ratio  tends 
to  be  largest  for  earthquakes,  somewhat  smaller  for  rockbursts,  and  much 
smaller  for  blasts;  however,  there  was  some  evidence  of  frequency 
dependence.  The  individual  phase  spectra  show  a  tendency  for  the  blast  P  to 
be  richest  in  high  frequencies  and  for  the  earthquake  P  to  be  enriched  at  high 
frequencies  relative  to  the  rockburst.  The  available  far-regional  and  teleseismic 
signals  from  larger  Central  European  and  South  African  events  showed  a 
similar  tendency  for  P  to  be  enriched  at  high-frequencies  for  earthquakes 
relative  to  rockbursts.  This  trend  was  strongest  for  the  largest  Central  European 
events.  Evidence  of  differences  in  spectral  behavior  for  some  far-regional 
events  and  effects  of  differences  in  propagation  path  at  regional  and  far- 
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regional  stations  will  require  additional  investigation  before  drawing 
conclusions  about  source  variability  between  individual  rockbursts. 

With  regard  to  other  techniques  for  discriminating  rockbursts,  rockbursts 
are  like  explosions  in  that  they  occur  at  comparable  focal  depths;  so  depth- 
based  discrimination  techniques  would  not  be  expected  to  provide  any 
distinction.  However,  some  characteristics  indicate  that  rockbursts  should 
behave  as  earthquakes  and,  therefore,  might  be  distinguished  from  nuclear 
explosion  tests  using  the  same  techniques  used  to  separate  earthquakes.  In 
particular,  rockbursts  which  correspond  to  release  of  tectonic  stress  on  pre¬ 
existing  zones  of  weakness,  or  faults,  near  the  excavation  should  behave  like 
earthquakes.  In  fact,  we  found  examples  in  our  comparisons  of  regional  signals 
for  South  African  events  that  S/P  ratios  in  the  time  domain  intermingled  for 
rockbursts  and  earthquakes. 

However,  some  rockbursts  are  more  complex  and  may  include  non¬ 
earthquake  processes  such  as  rock  fall  or  pillar  failure.  We  found  examples  in 
both  the  Central  European  and  South  African  databases  where  the  teleseismic, 
and  also  regional,  P  signals  were  at  least  as  complex  for  rockbursts  as  for 
earthquakes  from  the  same  region.  Such  complexity  in  rockburst  signals  may 
be  useful  as  a  discriminant  and  should  be  investigated  further.  The  complexity 
seen  in  the  rockburst  source  along  with  the  emergent  character  found  in  many 
regional  P  signals  from  rockbursts  tends  to  mitigate  against  the  value  of  sense 
of  first  motion  as  a  discriminant.  In  South  Africa  we  found  examples  where  the 
sense  of  first  motion  appeared  to  change  on  regional  records  from  rockbursts  in 
the  same  source  area.  Consequently,  some  rockburst  events  are  likely  to 
produce  anomalous  first  motion  patterns. 
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One  other  discrimination  technique  which  appears  to  be  problematic  for 
rockbursts  is  Ms  versus  mb-  In  the  first  place,  many  of  the  events  are  small  and 
would  require  some  kind  of  regional  measurement  of  the  magnitudes. 
Furthermore,  the  evidence  we  analyzed  for  even  large  rockbursts  indicated  a 
tendency  for  rockbursts  to  generate  relatively  weak  Ms.  This  explosion-like 
behavior  for  the  Ms/mb  discriminant  needs  further  study  including  more  stations 
and  other  source  areas. 

Finally,  with  regard  to  the  potential  use  of  rockbursts  to  conceal 
clandestine  nuclear  testing,  there  is  evidence  in  the  mining  literature  that 
rockburst  activity  can  be  controlled  by  engineering  techniques.  Application  of 
finite  element  procedures  to  specific  mining  areas  enables  the  engineer  to 
accurately  predict  stress  conditions  in  the  surrounding  rock.  Case  histories 
indicate  that  addition  or  removal  of  supports  can  be  used  to  reduce  stress  level 
below  some  threshold  or  increase  it  to  induce  a  rockburst  in  a  particular 
magnitude  range.  A  viable  evasion  scenario  appears  to  be  that  a  small  or 
decoupled  nuclear  explosion  test  might  be  concealed  in  a  triggered  rockburst. 
Additional  analyses  are  needed  to  determine  the  maximum  size  of  events  which 
can  be  controlled  in  this  manner.  If  large  rockbursts  can  be  controlled,  a 
simultaneously  detonated  nearby  nuclear  explosion  would  be  difficult  to  detect 
and  identify. 
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